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Abstract 


The  solubility  of  hydrogen  in  heavy  hydrocarbons  and 
coal-heavy  oil  slurries  was  measured  using  a  batch 
autoclave.  The  solvents  ranged  from  a  hydrocracked  gas  oil 
product  of  molecular  weight  250  to  heavy  bitumen  residues  of 
molecular  weights  above  800.  The  experiments  were  conducted 
mainly  at  200°  and  300°C,  and  pressures  up  to  24.8  MPa. 
Limited  carbon  dioxide  and  hydrogen  sulphide  solubility 

measurements  were  also  made. 

The  prediction  and  correlation  of  the  solubility  data 
by  the  Peng-Robinson  equation  and  the  Grayson-St reed  method 
were  attempted.  The  best  result  for  the  correlation  of 
hydrogen  solubility  data  was  obtained  with  the  use  of  either 
a  modified  Peng-Robinson  or  Soave-Redl ich-Kwong  equation  of 
state . 
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1  .  Introduction 


Solubility  data  for  gases  in  hydrocarbon  solvents  are 
useful  in  the  petroleum  industry.  The  effect  of  temperature 
on  gas  solubilities  is  of  practical  importance  because  of 
the  wide  range  of  temperatures  encountered  in  chemical 
processes.  The  measurement  and  correlation  of  phase 
equilibrium  data  provide  important  information  needed  in  the 
design  of  chemical  processing  equipment. 

Recent  interest  in  the  recovery  and  processing  of  shale 
oil,  coal  liquids  and  oil-sands  bitumens  has  created  a  need 
for  solubility  data  for  hydrogen,  carbon  dioxide  and 
hydrogen  sulphide,  at  elevated  temperatures  and  pressures, 
in  heavy  hydrocarbons  quite  different  from  the  normal 
petroleum  mixtures.  These  systems  contain  mono-nuclear  and 
poly-nuclear  aromatics,  naphthenes,  high  boiling  paraffins 
and  ring  compounds  containing  N,  S  and  0.  As  well,  the  need 
to  increase  recovery  of  distillates  from  petroleum  residues 
has  demanded  accurate  vapor-liquid  equilibrium  data  for 
hydrogen  in  heavy  petroleum  mixtures. 

In  hydrotreating  applications  hydrogen  coexists  with 
carbon  dioxide,  hydrogen  sulphide  and  light  hydrocarbon 
gases  in  hydrocarbon  liquids.  The  design  of  reformers, 
hydrocrackers,  hydrotreaters  and  scrubbing  units  requires 
accurate  knowledge  of  solubilities  of  these  gases. 

The  Canada  Centre  for  Mineral  and  Energy  Technology 
(CANMET)  of  the  Department  of  Energy,  Mines  and  Resources 
has  developed  a  hydrocracking  process,  operable  at  high 


1 


2 


pressure  and  producing  low  pitch  and  high  distillate  yields 
from  bitumen  fractions.  Figure  1  presents  a  simplified  flow 
scheme  of  a  one-barrel  per  day  CANMET  hydrocracking  pilot 


plant  simulating  a  commercial  process.  A  mixture  of  hydrogen 
and  topped  bitumen  is  preheated  and  pumped  into  a  4.5  litre 
reactor  vessel  where  the  bitumen  is  cracked  and  hydrogenated 
to  form  low  molecular  weight  distillate  and  residual  pitch. 
The  reactor  effluent  is  cooled  in  two  stages  giving  a  heavy 
oil  and  a  light  oil.  The  hydrogen-rich  gas  stream  is 
scrubbed  to  remove  gaseous  hydrocarbons  and  hydrogen 
sulphide  ,  mixed  with  make-up  hydrogen  and  recycled  back  to 
the  reactor.  Because  of  the  high  sulphur  content  in  bitumen, 
substantial  hydrogen  sulphide  is  formed  during 
hydrocracking.  The  process  is  capable  of  handling  a  wide 
variety  of  feedstocks.  The  pitch  conversion  can  be  altered, 

as  desired. 

Hydrogen  exists  in  equilibrium  with  the  feed  and 
products  at  the  high  temperature  required  in  the  CANMET 
hydrocracking  process,  together  with  hydrogen  sulphide 
formed  by  hydrodesulphur i sat  ion .  Consequently,  solubility 
data  for  hydrogen  and  hydrogen  sulphide  in  bitumen 
feedstocks  and  cracked  products  are  useful  for  the  design 
and  analysis  of  commercial  bitumen-upgrading  processes.  As 
well,  these  data  are  required  in  the  kinetic  modelling  of 
the  complex  reactions  taking  place  during  the  hydrocracking 
and  hydrodesulphur i sat i on .  The  present  work  was  undertaken 
to  provide  data  at  the  temperatures  and  pressures 
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Schematic  Diagram  of  Canmet  Pilot  Plant 
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approaching  those  used  in  the  CANMET  hydrocracking  process. 
The  solvents  used  were  actual  feedstocks  and  hydrocracked 
products  from  the  pilot  plant. 

Coal  upgrading  can  be  done  in  a  slurry  medium;  the 
solvent  used  can  either  donate  hydrogen  or  act  as  a  hydrogen 
carrier.  No  data  on  hydrogen  solubility  in  coal-oil  slurries 
are  available.  These  studies  were  extended  to  obtain 
hydrogen  solubility  in  selected  coal-bitumen  slurries.  The 
data  will  be  helpful  in  the  analysis  of  kinetics  of  coal 
liquefaction . 

As  a  result  of  the  present  shortage  of  crude  oil,  there 
is  increased  interest  in  the  use  of  carbon  dioxide  as  a 
tertiary  oil  recovery  agent,  since  it  promotes  swelling  of 
the  heavy  oils  and  reduces  their  viscosity  greatly.  The 
knowledge  of  carbon  dioxide  solubility  in  bitumens  at 
reservoir  temperatures,  as  well  as  at  temperatures  required 
to  recover  the  solute  gas,  is  desirable. 

In  the  modelling  of  hydrocracking  and  related 
hydrotreatment  processes,  thermodynamic  correlations  are 
used  to  calculate  the  vapor  and  liquid  phase  compositions 
and  their  properties.  The  ’dual'  model  method  of  Chao-Seader 
is  a  technique  widely  applied  in  the  petroleum  industry.  The 
Soave  modification  of  the  Redl ich-Kwong  equation  of  state 
and  the  Peng-Robinson  equation  of  state  are  more  recent 
methods,  and  are  popular  due  to  their  simplicity  and 
accuracy.  The  reliability  of  these  methods  to  predict  phase 
equilibrium  was  tested  using  experimental  data  obtained  here 


and  from  the  literature. 
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2.  Literature  Survey 


2 . 1  Theory 

For  most  gases  at  ordinary  temperatures  the  solubility 
in  a  liquid  solvent  decreases  with  rising  temperature.  This 
is  expected  because,  as  temperature  rises,  the  kinetic 
energy  of  the  gas  increases  and  thus  there  is  a  decreasing 
tendency  to  condense  into  a  liquid  phase.  However,  the 
solubilities  of  some  gases  increase  with  temperature  while 
others  first  decrease  and  then  increase.  It  appears  likely 
that  the  solubilities  of  all  gases  first  decrease,  go 
through  a  minimum,  and  then  increase  if  taken  over  a  wide 
enough  temperature.  The  increase  of  solubility  of  some  gases 
with  temperature  is  probably  related  to  the  decrease  in 
solvent  density  as  temperature  rises  (45,69). 

The  temperature  derivative  of  the  solubility, 
calculated  from  the  Gibbs- Helmhol tz  equation,  is  directly 
related  to  the  partial  molar  entropy  of  the  gaseous  solute 
in  the  liquid  phase.  Therefore,  if  something  can  be  said 
about  the  entropy  change  of  the  solution,  insight  can  be 
gained  on  the  effect  of  temperature  on  solubility.  Consider 
the  relatively  simple  case  where  the  solvent  is  essentially 
non-volatile  and  where  the  solubility  is  sufficiently  small 
to  make  the  activity  coefficient  of  the  solute  independent 
of  the  mole  fraction.  With  these  restrictions  it  can  be 

shown  that 
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3  In  x2  As 2 


(2.1) 


3  In  T 


R 


where  x2  is  the  mole  fraction  of  gaseous  solute  at 
saturation  and 


(2.2) 


As2  =  s2 (1)  -  s2 (v) 


First  we  consider  equation  (2.1);  if  the  partial  molar 
entropy  of  change  of  the  solute  is  positive,  then  the 
solubility  increases  with  rising  temperature;  otherwise  it 
falls.  To  understand  the  significance  of  the  entropy  change 
it  is  convenient  to  divide  it  into  two  parts; 


(2.3) 


As 2  =  (s2(l)  -  s2(v))  +  ( s 2 ( 1 )  ~  s2(l)) 


where  s2(l)  is  the  entropy  of  the  (hypothetical)  pure  liquid 
at  the  temperature  of  the  solution.  The  first  term  on  the 
right-hand  side  of  equation  (2.3)  is  the  entropy  of 
condensation  of  the  pure  gas  and  in  genera-  this  is  expec  _ ec 
to  be  negative  since  the  entropy  of  a  liquid  is  lower  than 
that  of  a  saturated  gas  at  the  same  temperature.  The  secona 
term  is  the  partial  molar  entropy  of  solution  of  the 
condensed  solute  and,  assuming  idea-  entropy  c _  mixing  -3- 
the  two  liquids,  it  is 


s2  ( 1 )  -  s 2 ( 1 ) 
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(2.4) 


i 


8 


Since  x2<1  the  second  term  in  Eq  (2.3)  is  positive  and 
the  smaller  the  solubility,  the  larger  this  term  becomes.  It 
therefore  follows  that  As2  should  be  positive  for  those 
gases  which  have  very  small  solubilities  and  negative  for 
others.  Gases  which  are  sparingly  soluble  show  positive 
temperature  coefficients  of  solubility,  whereas  gases  which 
are  readily  soluble  show  negative  temperature  coefficients. 
This  is  in  fact  observed. 


2.2  Experimental  Techniques 

Solubility  determinations  consist  essentially  in 
bringing  the  solvent  into  equilibrium  with  the  gas  at  a 
known  total  pressure,  taking  a  sample  of  the  liquid  phase 
and  determining  its  composition.  Individual  methods  differ 
in  the  manner  in  which  the  attainment  of  equilibrium  is 
accelerated  ( stirring , shaking  circulating  the  liquid  or 
bubbling  the  gas)  and  as  to  whether  sampling  is  continuous 
or  intermittent.  Chappelow  and  Prausnitz  (5) ,  and  Cukor  and 
Prausnitz  (8)  investigated  the  solubility  of  hydrogen  at  low 
pressure  in  bicyclohexyl,  dipheny lmethane ,  squalane, 
n-hexadecane  and  octamethyl-cyclotetras i loxane .  The 
experimental  apparatus  consisted  of  5  sections.  (1)  a 
degassing  flask,  (2)  an  equilibrium  cell,  (3)  a  reservoir  of 
degassed  solvent,  (4)  a  precision  gas  burette,  and  (5)  a 
section  for  pressure  measurement.  Solvent  was  continuously 
circulated  through  a  vapor  space  containing  a  measured 
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quantity  of  solute.  The  solubility  was  determined  by  making 
a  material  balance  on  the  gaseous  solute.  The  precision  of 
their  results  was  about  1%. 

The  solubilities  of  hydrogen  in  isobutane, 

2 , 2 , 4-tr imethylpentane ,  and  in  a  mixture  of  isomeric 
dodecanes  were  investigated  by  Dean  and  Tooke  ( 1 1 ) .Hydrogen 
and  the  liquids  were  placed  in  an  equilibrium  cell  mounted 
in  an  oil-filled  bath  in  which  it  was  rocked  at  30  cycles  a 
minute.  At  the  upper  end  of  the  cell  one  valve  was  attached 
with  an  induction  tube  for  withdrawing  liquid  phase  samples. 
Another  valve  was  provided  for  taking  vapor  phase  samples. 
Solubilities  in  isobutane  were  determined  for  temperatures 
from  100°  to  250°F  and  pressures  from  500  psia  to  3000  psia. 
The  solubility  of  hydrogen  in  2 , 2 , 4-tr imethylpentane  was 
measured  for  temperatures  from  100°  to  302. 5°F  and  in  the 
dodecanes  for  200°  and  300°F  with  pressures  ranging  from  500 
to  5000  psia.  The  liquid  and  vapor  samples  were  analyzed  by 
repeatedly  freezing  and  thawing  the  samples,  resulting  in 
the  physical  separation  of  the  hydrogen  from  the 
hydrocarbon.  An  Orsat  apparatus  was  also  used  to  analyze  the 
hydrogen-isobutane  systems.  They  showed  that  hydrogen 
solubility  increased  with  temperature  and  pressure  ana  that 
solubility  decreased  as  solvent  molecular  weight  increased. 
They  also  concluded  that  hydrogen  is  more  soluble  in 
paraffins  than  in  aromatics  of  similar  molecular  weight. 

Eakin  and  DeVaney  (12)  investigated  the  vapor-liquid 
equilibrium  for  hydrogen  sulphide-n-nonane ,  hydrogen 
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sulphide- i sopropylcyc lohexane ,  and  hydrogen  sulphide- 
mesitylene  at  100°  to  400°F.  Ternary  mixtures  of  hydrogen 
sulphide-hydrogen-n-nonane ,  hydrogen  sulphide-hydrogen- 
i sopropylcyc lohexane ,  and  hydrogen  sulphide-hydrogen- 
mesitylene  at  up  to  2000  psia  at  these  same  temperatures 
were  also  studied.  The  apparatus  consisted  of  a  rocking  type 
equilibrium  cell.  All  samples  were  analyzed  on  a  Beckman 
GC-4  chromatograph.  The  column  was  a  16  foot  long,  1/8  inch 
OD  tube  packed  with  80-100  mesh  Porapak  P.  The  carrier  gas 
was  argon. 

Frolich  et  al.  (15)  determined  gas  solubility  by 
saturating  a  liquid  with  gas  at  suitable  intervals  of 
pressure  and  measuring  the  ratio  of  solute  to  solvent  in  a 
sample  drawn  off  for  analysis.  To  this  end  a  small  amount  of 
solvent  was  introduced  into  an  evacuated  steel  cylinder  of  2 
litre  capacity.  The  gas  was  then  forced  in  at  the  highest 
pressure  available  and  the  cylinder  agitated  in  a  water  bath 
maintained  at  25°C.  The  liquid  phase  analysis  was  done  by 
flashing  a  sample  of  the  liquid  and  collecting  the  vapor 
over  mercury  in  a  burette.  The  total  error  of  this  method  is 
such  that  the  results  should  not  be  considered  accurate  to 
better  than  ±5%.  The  liquids  studied  included:  gas  oil, 
heavy  naphtha,  cyclohexane,  octane,  and  hexane  at  25°C  and 

pressures  up  to  200  atmospheres. 

Grove  et  al.  (19)  measured  the  solubility  of  hydrogen 
in  terphenyl  at  high  temperatures.  Their  apparatus  was  made 
of  two  cylindrical  stainless  steel  vessels  isolated  from  one 
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another  by  means  of  a  stainless  steel  bellows.  A  pressure 
difference  between  the  two  vessels  caused  contraction  or 
expansion  of  the  bellows.  The  assembly  was  agitated'  by  a 
motor  at  about  320  cycles/min  to  establish  equilibrium.  A 
known  weight  of  terphenyl  was  placed  in  one  vessel,  degassed 
and  pressured  with  hydrogen,  the  other  vessel  being 
pressured  with  nitrogen  simultaneously  so  that  no  large 
stresses  were  placed  on  the  bellows.  Hydrogen  solubilities 
were  calculated  by  measuring  the  decrease  in  the  number  of 
moles  of  gas  in  the  vapor  space  over  the  solvent  in  passing 
from  the  initial  equilibrium  condition  when  the  solvent  was 
solid,  to  the  temperature  of  measurement. 

Ipatieff  et  al.  (21)  studied  the  solubility  of  hydrogen 
in  benzene,  toluene,  xylene,  mesitylene,  gasoline,  kerosene, 
cylinder  oil  and  paraffin  mazout .  The  liquid  under 
investigation  was  placed  in  a  high  pressure  bomb.  The 
necessary  pressure  was  applied  through  an  upper  valve  and 
the  liquid  was  thoroughly  mixed  with  the  gas.  Afterwards  a 
sample  was  taken  into  a  gas  burette  through  a  special  valve. 
The  volumes  of  the  liquid  and  the  evolved  gas  were  read 
directly  from  the  burette.  Solubilities  were  reported  for 
temperatures  of  25,  40  and  100°C  and  for  pressures  up  to  319 
atmospheres  in  benzene,  toluene,  xylene,  and  mesitylene,  and 
at  25°  to  300°C  (to  a  maximum  pressure  of  300  atmospheres) 
in  gasoline  fraction,  kerosene,  cylinder  oil  and  paraffin 
mazout  oil.  They  concluded  that  (i)  solubilities  increased 
as  temperature  increased  (ii)  solubilities 
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linearly  with  pressures  to  300  atmospheres  and  (iii) 
solubilities  decreased  with  the  complexity  of  the 
composition  of  the  molecule,  being  the  greatest  for  benzene 
and  gasoline,  and  the  smallest  for  mesitylene  and  mazout. 

Lachowicz  et  al.  (26)  investigated  the  solubility  of 
hydrogen  in  n-heptane  and  n-octane  at  25°,  37.5°  and  50°C  at 
pressures  between  50  and  300  atmospheres.  Their  apparatus 
consisted  of  (i)  a  mercury  piston  compressor,  (ii)  an 
equilibrium  cell,  and  (iii)  an  apparatus  for  expanding  the 
sample  to  1  atm.  and  analyzing  it.  The  vessel  was  rocked  in 
a  thermostat  to  attain  equilibrium.  During  the  withdrawal  of 
the  sample,  the  pressure  in  the  vessel  was  maintained  within 
±0.5  atm.  by  admitting  compressed  gas  from  a  gas  storage 
vessel.  The  liquid  sample  was  expanded  into  a  cooled  trap 
and  frozen.  The  evolved  gas  was  collected  in  a  burette  and 
its  volume  measured  at  room  temperature  and  pressure.  The 
liquid  sample  was  weighed  after  bringing  it  to  room 
temperature . 

Nichols  et  al.  (39)  measured  the  specific  volumes  of 
four  mixtures  of  hydrogen-/?- hexane  at  eight  temperatures 
between  40  and  460°F,  for  pressures  up  to  10000  psi.  A 
mixture  sample  of  known  composition  and  weight  was  confined 
in  a  stainless  steel  vessel  over  mercury.  The  effective 
volume  of  this  chamber  was  varied  by  the  introduction  and 
withdrawal  of  mercury.  Mechanical  agitation  was  provided  to 
hasten  the  attainment  of  physical  equilibrium  within  and 
between  the  phases.  The  cell  was  immersed  in  an  agitarec 
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liquid  bath.  Sample  of  the  gas  phase  was  removed  under 
i sobar ic- i sothermal  conditions  and  passed  through  a  series 
of  weighing  bombs  maintained  at  liquid  nitrogen  temperature. 
The  total  quantity  of  hydrogen  in  the  gas  phase  was  then 
determined  volumetr ically .  No  liquid  compositions  were 
reported . 

Prather  et  al.  (44)  used  an  Autoclave  Engineers  316 
stainless  steel  Magnedrive  autoclave  with  a  volume  of  one  US 
gallon  to  study  the  solubility  of  hydrogen  in  recycle  oil 
and  creosote  oil,  and  in  coal-creosote  oil  slurry.  Hydrogen 
analysis  was  carried  out  on  a  Varian  model  920  Aerograph 
using  a  15  ft  column  packed  with  75%  molecular  sieve  13  X 
and  25%  molecular  sieve  5A.  The  oil  was  charged  into  the 
autoclave  and  the  system  was  evacuated.  The  autoclave  was 
then  brought  to  the  desired  temperature  and  hydrogen  was 
added.  Stirring  was  carried  on  at  2000  rpm.  Samples  of  the 
vapor  were  withdrawn  from  the  top  of  the  autoclave  and 
analyzed  by  GC.  Liquid  samples  were  withdrawn  from  the 
autoclave  into  a  stainless  steel  bomb.  The  bomb  was  then 
fitted  to  an  evacuated  glass  rack,  in  which  the  volume  of 
the  dissolved  gases  was  measured.  The  gas  was  withdrawn  by 
action  of  a  Toepler  pump,  and  after  volume  measurement  it 
was  passed  to  a  GC  for  analysis.  Data  were  provided  at 
temperatures  up  to  400°C  and  for  385  psia  to  3000  psia 
hydrogen  pressure.  These  authors  concluded  that  hydrogen 
solubility  in  a  3:1  solvent:coal  mixture  was  similar  to  that 
in  the  pure  solvent  at  temperatures  high  enough  for  coal  to 
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be  thermally  decomposed. 

Prausnitz  and  Benson  (46)  studied  the  vapor  phase 
solubility  of  iC8H18,  C6H5CH3,  and  C10H22  in  compressed 
hydrogen  with  a  flow  system.  The  gas  was  dried  with  silica 
gel  before  entering  the  equilibrium  cells.  Two  equilibrium 
cells  connected  in  series  were  used  to  assure  saturation  of 
the  gas.  The  liquid  was  contained  within  the  cell,  in  pyrex 
liners,  and  the  gas  was  very  finely  dispersed  from  a  fritted 
glass  sparger  (pore  size  14  microns)  at  the  bottom  of  the 
cells.  The  equilibrium  cells  were  immersed  in  a  constant 
-temperature  bath,  which  was  controlled  to  within  ±0.05°C. 
After  leaving  the  second  cell  the  vapor  was  led  to  a  needle 
expansion  valve.  After  expansion  the  vapor  mixture  entered 
the  condenser.  It  was  assumed  that  the  vapor  leaving  the 
condenser  was  in  equilibrium  with  the  condensate.  The  flow 
of  gas  was  measured  with  a  wet-test  meter.  From  the  amount 
of  condensate  collected  and  hydrogen  passed  through  the 
wet-test  meter,  the  vapor  phase  composition  was  determined. 

Simnick  et  al.  (61,62,63,64),  Sebastian  et  al. 
(53,55,56),  Lin  (29),  and  Yao  et  al.  (72)  studied  the 
gas-liquid  equlibrium  for  binary  mixtures  of  hydrogen  in 
/77-cresol,  /77-xylene,  tetralin,  dipheny  lmethane ,  bicyclohexyl, 
toluene,  1 -methy lnaphthalene ,  thianaphthene  and  quinoline.  A 
flow  apparatus  was  used  to  reduce  the  residence  time  of  the 
fluids  in  the  apparatus.  The  liquid  was  pumped  in  at  a  rate 
of  500-2000  cc/h  and  the  hydrogen  gas  was  supplied  at  a  rate 
of  about  5500  cc/min.  The  two  streams  were  mixed,  and  heated 


■ 


15 


in  a  tubular  heater  to  the  desired  temperature  and  then 
passed  into  an  isothermal  cell.  The  cell  effluent  was 
throttled  and  cooled  to  ambient  conditions.  The  liquid  was 
trapped  and  later  weighed  with  an  analytical  balance  and  the 
quantity  of  hydrogen  gas  determined  volumetr ically .  The  work 
was  later  extended  to  ternary  systems  of  hydrogen 
-tetral in-diphenylmethane  and  hydrogen-tet ral in-m-xylene 
mixtures . 

Sagara  et  al.  (51)  determined  the  vapor-liquid 
equilibrium  relations  of  the  hydrogen-methane,  hydrogen- 
ethylene,  hydrogen-ethane,  hydrogen-methane-ethylene  and 
hydrogen-ethane-ethylene  systems  by  a  static  method  at 
several  temperatures  from  -25°  to  -170°C  and  at  pressures  up 
to  100  atmospheres.  The  vapor  and  liquid  samples  were 
collected  in  chambers  and  analyzed  by  gas  chromatography. 

The  work  which  reached  the  highest  temperature  was  by 
Grayson  and  Streed  (18)  who  presented  a  correlation  of  K 
values  of  hydrogen  in  heavy  oils  at  320°  to  480°C  and 
pressures  to  200  atmospheres.  No  experimental  data  were 

reported . 

King  et  al.  (24,25)  measured  the  solubilities  of  carbon 
dioxide  and  hydrogen  sulphide  as  functions  of  temperature  in 
the  r?-alkanes  hexane  to  hexadecane  at  normal  pressure,  and 
correlated  the  solubilities  at  25°C.  A  'falling  film'  type 
of  apparatus  was  used.  Each  solvent  was  degassed  before  use 
by  spraying  into  a  continuously  evacuated  chamber.  It  was 
then  displaced  from  this  chamber  by  injecting  mercury  and 
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flowed  drop  by  drop  over  a  lip  where  it  passed  as  a  film 
down  the  inside  wall  of  a  spiral,  and  finally  accumulated  in 
a  graduated  pipette.  A  magnetically  operated  stirrer 
agitated  the  gas-liquid  interface  in  the  pipette.  A  known 
volume  of  gas  was  stored  in  a  burette  and  amount  dissolved 
in  a  known  quantity  of  solvent  was  determined  by  noting  the 
volume  of  gas  expelled  from  the  gas  burette. 

Low  pressure  solubilities  for  hydrogen  sulphide  in 
/7-hexadecane ,  dipheny  lmethane ,  bicyclohexyl,  and 
1 -methylnaphthalene  have  been  studied  by  Tremper  and 
Prausnitz  (69).  The  experimental  apparatus  and  procedure  was 
essentially  the  same  as  those  described  by  Cukor  and 
Prausnitz  ( 9  )  . 

Yarborough  (73)  investigated  the  vapor-liquid 
equilibrium  for  hydrogen  sulphide  and  n~ decane.  An 
equilibrium  cell,  sampling  valves,  and  2  different 
chromatographs  were  used  for  analysis  during  this  study. 

Jacoby  and  Rzasa  (22)  obtained  experimental 
vaporization  K-values  for  hydrogen  sulphide  and  carbon 
dioxide  in  two  natural  gas-absorber  oil  mixtures  and  in  two 
natural  gas-crude  oil  mixtures  at  temperatures  of  100,  150 
and  2 0 0  0 F  and  various  pressures  in  the  range  200  to  5000 
psia.  The  pressure  cell  consisted  of  a  hollow  stainless 
steel  cylinder,  placed  horizontally,  and  fitted  with  a 
piston.  The  liquid  phase  was  circulated  with  a  magnetic 
pump.  The  vapor  samples  were  analyzed  directly  with  a  mass 
spectograph.  The  liquid  samples  were  flashed  and  the  light 
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gases  were  analyzed  by  the  mass  spectrograph.  The  remainder 
of  the  samples  were  fractionated  and  C7+  fractions  were 
character i sed . 

Simon  and  Graue  (66)  studied  the  solubility  of  carbon 
dioxide  in  nine  different  oils  at  conditions  covering  a 
range  of  100  to  250  °F  and  pressures  up  to  2300  psia.  Known 
amounts  of  pure  carbon  dioxide  and  crude  oil  were  combined 
in  a  visual  cell  and  bubble  points  of  the  mixture  were 
measured.  Correlations  for  predicting  the  solubility, 
swelling,  and  viscosity  behaviour  of  these  systems  were 
obtained  in  graphical  form. 

From  the  experimental  methods  described  by  the  various 
authors,  it  was  concluded  that  a  batch  autoclave  is  suitable 
for  the  measurement  of  hydrogen  solubility  in  heavy 
bituminous  materials  and  coal-bitumen  slurries.  This  type  of 
equipment  is  available  in  different  sizes,  and  operable  at 
high  temperatures  and  pressures.  A  GC ,  with  either  N2  or  Ar 
as  carrier  gas,  is  adequate  for  both  the  liquid  and  vapor 
phase  analyses. 


2.3  Correlation 

Early  attempts  to  describe  the  volumetric  properties  of 
gases  led  to  the  ideal  gas  law  and  the  van  der  Waals 
equation  for  real  gases.  Since  then  systematic  efforts  have 
been  made  to  describe  the  equilibrium  properties  of  real 
gases.  A  review  of  the  various  equations  oi  state  that  have 
been  widely  used  is  presented  by  Tsonopoulos  and  Prausnitz 
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(70)  . 


Equations  of  state  serve  four  purposes  :  (i)  prediction 
of  vapo'r-1  iquid  equilibria  of  mixtures  at  wide  ranges  of 
temperatures  and  pressures,  (ii)  prediction  of  gas  phase 
properties  of  pure  fluids  and  their  mixtures,  (iii) 
prediction  of  liquid  phase  properties  of  pure  fluids  and 
their  mixtures,  (iv)  interpolation  of  PVT  data,  and 
differentiation  and  integration  of  PVT  data  for  estimation 
of  derived  properties  (e.g.  fugacity). 

A  single  equation  of  state  does  not  satisfactorily 
achieve  all  these  purposes.  The  virial  equation  of  state 
expresses  the  compressibility  factor,  z,  of  a  gas  by  a  power 
series  in  the  density,  and  is  useful  only  for  modest 
densities  (approximately  3/4  of  the  critical  density)  when 
virial  coefficients  beyond  the  third  are  ignored.  The  virial 
equation  is  attractive  because  of  its  theoretical  basis  and 
of  its  relation  to  intermolecular  forces.  Its  usefulness  is 


limited  by  our  incomplete  knowledge  of  virial  coef f iceients . 

Benedict,  Webb  and  Rubin  (2)  presented  an  equation  of 
state  with  eight  adjustable  parameters.  To  obtain  eight 
meaningful  constants,  extensive  experimental  data  are 
required.  Even  if  such  data  are  available,  the  constants  are 
not  unique.  The  equation  represents  satisfactorily 
volumetric  properties  of  non-polar  gases  and  liquids  to 
densities  up  to  about  1.8  times  the  critical. 

The  most  successful  variation  on  the  van  der  Waals 
equation  was  that  proposed  by  Redlich  and  Kwong  (49),  which 
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uses  only  two  adjustable  parameters.  The  RK  equation,  while 
admittedly  less  accurate  than  the  BWR  equation  when  applied 
to  pure  substances,  is  attractive  for  mixtures  since  it  does 
not  rely  on  extensive  experimental  data.  The  equation  has 
been  applied  extensively  to  vapor  mixtures  and  their 
properties.  A  notable  application  is  in  the  Chao-Seader  (4) 
correlation  which  applies  a  ’dual'  model  approach  to 
calculate  the  vapor-liquid  equilibria.  The  vapor-phase 
fugacity  coefficient  is  calculated  using  the  RK  equation  of 
state . 

Numerous  variations  on  the  RK  equation  have  been 
proposed.  Chueh  and  Prausnitz  (6)  applied  the  RK  equation 
assuming  a  geometric  mean  mixing  rule  for  the  critical 
temperature  of  the  binary  formed  by  i  and  j,  and  introducing 
the  interaction  parameter,  kjj,  to  correct  for  the  departure 
from  the  geometric  mean. 

Zudkevitch  and  Joffe  (74)  also  modified  the  RK 
equation,  utilising  Chueh  and  Prausnitz' s  ideas,  and 
introduced  the  binary  interaction  parameter  to  correct  for 
the  deviation  from  the  geometric  mean  in  the  cross  parameter 
a  j  j .  They  applied  the  RK  equation  to  both  the  vapor  and  the 

liquid  phases. 

In  the  Chao-Seader  correlation  the  upper  limit  for  the 
liquid  fugacity  coefficient  for  hydrogen  is  260  C.  Grayson 
and  Streed  (18)  used  Chao  and  Seader's  experimental  data,  as 
well  as  their  own,  to  modify  the  Chao-Seader  correlation  for 
the  K-values  in  hydrogen-gas  oil  systems.  The  correlation  is 
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applicable  for  temperatures  up  to  480°C. 

Recently,  more  emphasis  has  been  placed  on  applying  a 
single  equation  of  state  for  the  calculation  of  vapor-liquid 
equilibria,  since  the  same  equation  can  then  be  used  for  the 
fugacity  coefficient  in  both  the  phases.  Soave  (67)  modified 
the  attraction  parameter  term  in  the  RK  equation,  making  it 
temperature  dependent,  as  well  as  a  function  of  the  acentric 
factor.  The  SRK  equation  of  state  gained  acceptance  in  the 
hydrocarbon  processing  industry  because  of  its  relative 
simplicity  and  its  ability  to  generate  reasonably  accurate 
K-values  in  vapor-liquid  equilibrium  calculations.  However, 
both  the  SRK  and  RK  equations  fail  to  generate  satisfactory 
density  values  for  the  liquid.  Graboski  and  Daubert  (17) 
recommended  improvements  to  the  SRK  method  to  handle 
hydrogen-containing  systems. 

Peng  and  Robinson  (41)  presented  a  simple  equation  of 
state  that  gave  improved  liquid  density  values  as  well  as 
accurate  vapor  pressures  and  K— values.  The  PR  method  assumed 
the  geometric  mean  mixing  rule  in  the  energy  parameter,  a  , 
for  the  binary.  An  interaction  coefficient  was  then  applied 
to  correct  for  the  departure  from  the  geometric  mean.  The 
one-fluid  approximation  was  used  to  define  the  covolume 
factor,  b,  for  the  mixture. 

El-Twaty  and  Prausnitz  (14)  obtained  an  improved 
correlation  of  K-values  for  hydrogen-heavy  hydrocarbon 
mixtures  by  modifying  the  composition  dependence  of  the 
covolume  factor  b  in  the  SRK  equation.  The  modification 
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introduced  a  binary  interaction  parameter  EHj  that  can  be 
adjusted  to  fit  binary  data. 

Recently,  Sebastian  et  al.  (54)  developed  a  correlation 
in  which  the  fugacity  of  dissolved  hydrogen  at  zero  pressure 
was  correlated  as  a  function  of  temperature  and  solubility 
parameter.  The  high  pressure  fugacity  was  obtained  upon 
applying  a  Poynting  factor.  The  new  correlation  has  its 
basis  in  the  'regular’  solution  assumption,  and  the 
Chueh-Prausn i tz  modification  of  the  RK  equation  of  state. 

A  literature  search  has  shown  that  many  equations  of 
state  have  been  proposed  for  mixtures  containing  low 
molecular  weight  solvents.  Only  a  few  of  these  have  been 
applied  to  highly  asymmetric  mixtures  of  small  molecules, 
such  as  hydrogen,  and  large  molecules,  such  as  tetralin. 
This  work  considers  the  applicability  of  the  PR,  the 
modified  SRK  and  the  Grayson-St reed  methods  for  predicting 
solubi 1 i ty  of  hydrogen  in  high  boiling  point  hydrocarbon 
solvents  as  well  in  the  heavier  and  more  complex  feedstocks 
and  products  of  the  CANMET  bitumen  upgrading  process.  The 
ability  of  the  PR  equation  of  state  to  handle  hydrogen 
sulphide  and  carbon  dioxide  in  these  solvents  was  also 

investigated . 

2.3.1  Grayson-Streed 

The  Chao-Seader  correlation  is  limited  to  a  maximum 
temperature  of  260°C  for  hydrogen  systems.  New  experimental 
vapor-liquid  equilibria  data  for  high  temperature,  high 
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pressure  hydrogen-hydrocarbon  systems  were  used  to  extend 
the  temperature  range  of  the  Chao-Seader  K  correlation.  The 
experimental  results  covered  a  range  of  temperature  up  to 
480°C  and  a  range  of  pressure  up  to  3000  psia.  Modified 
coefficients  for  computing  the  liquid  fugacity  coefficient 
at  temperatures  up  to  425°C  for  hydrogen,  methane  and 
heavier  hydrocarbons  were  developed. 

In  the  Chao-Seader  correlation  K-values  are  calculated 
through  a  combination  of  three  factors: 


K  =  y  /  x  =  ^ !  y  /  <t> 


(2.5) 


The  quantity  vy  is  the  fugacity  coefficient  of  the  pure 
component  in  the  liquid  phase.  The  quantity  y  is  the 
activity  coefficient  of  this  component  in  the  liquid 
mixture.  The  quantity  <fi  is  the  fugacity  coefficient  of  the 
component  in  the  vapor  mixture. 

The  upper  temperature  limits  in  the  Chao-Seader 
correlation  for  the  liquid  fugacity  coefficient,  v\,  are: 

1.  Hydrogen  and  methane:  260°C 

2.  Other  hydrocarbons:  reduced  temperature  of  1.3. 

The  liquid  activity  coefficient  and  vapor  fugacity 
coefficient  correlations  remain  unchanged  from  those 
proposed  by  Chao-Seader.  Liquid  activity  coefficients  are 
calculated  from  Hildebrand's  equation,  assuming  regular 
liquid  solutions.  Vapor  fugacity  coefficients  are  calculated 
from  the  Redl ich-Kwong  equation  of  state. 
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The  Chao-Seader  correlation  for  liquid  fugacity 
coefficient  has  been  developed  within  the  framework  of  Curl 
and  Pitzer's  modified  form  of  the  principle  of  corresponding 
states  which  is  presented  in  equation  (2.6): 

log  Vy  =  log  v  2  +  gj  log  (2.6) 

The  first  term  on  the  right-hand  side  represents  the 
fugacity  coefficient  of  ’simple  fluids'.  The  second  term  is 
a  correction  accounting  for  the  departure  of  the  properties 
of  ’simple  fluids'. 

The  quantity  v2  is  dependent  only  on  reduced 
temperature  and  reduced  pressure  and  was  fitted  with  the 
following  function  by  Chao  and  Seader: 

log  v2  =  Ao  +  A^Tr  +  A2Tr  +  A3Tr  2  +  AuTr  3  +  (A5  +  A6Tr  + 
A 7 Tr  2  ) Pr  +  (A8  +  A9Tr  )Pr  2  -  log  Pr  (2.7) 

New  coefficients  for  equation  (2.7)  were  determined  to 
fit  Grayson  and  Streed's  experimental  high  temperature  data 
as  well  as  Chao  and  Seader's  low  temperature  data.  A 
separate  set  of  coefficients  are  required  for  hydrogen, 
since  the  typical  application  temperatures  in  the 
hydrotreating  processes  are  far  above  the  critical 
temperature  of  hydrogen.  These  are  listed  in  Grayson  and 
Streed  (18).  The  value  of  co ,  the  acentric  factor,  for 
hydrogen  is  then  taken  as  zero  in  equation  (2.6). 
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The  quantity  v3  is  also  dependent  only  on  reduced 
temperature  and  reduced  pressure  and  was  fitted  with  the 
following  function  by  Chao  and  Seader: 

log  =  -  4.23893  +  8.65808  Tr  -  1  .22060  /  Tr  - 

3  .  15224  Tr  3  -  0.025  (Pr  -  0.6)  (2.8) 

A  modification  to  the  limits  of  this  equation  is 
necessary  in  order  to  fit  the  high  temperature  data.  For 
values  of  Tr  above  1.0,  values  of  are  set  equal  to  the 
value  at  Tr  =  1.0. 


2.3.2  Soave-Redl ich-Kwong  Method 

Soave  (67)  originally  proposed  a  modification  of  the 
Redl ich-Kwong  equation  of  state  which  introduced  a  third 
parameter,  the  acentric  factor,  and  a  temperature  dependency 
into  the  cohesive  energy  term  to  account  for  the  effect  of 
non-spher ic i ty  on  fluid  P-V-T  properties.  The  forms  of  the 
equation  of  state  for  the  compressibility  factor  and 
fugacity  coefficient  are  given  by  equations  (2.9)  and 

(2.10): 


v  a  a 

v-b  RT ( v+b ) 


z 


(2.9) 


25 


In  0 


V 

bi 

z 

v+b 

In  ( - )  + 

(  ) 

-  -  I  a  j  j  x  j  In 

(  )  + 

v+b 

v-b 

RTb  j 

b 

ab  j 

-  {  In 

v+b 
(  ) 

b 

-  ( - )}  -  In  z 

(2.10) 

RTb 2 

b 

v+b 

For  any  component,  the  constants  a  and  b  are  determined 
from  the  universal  relationships: 


a  =  0.42747  R2  Tc  2  /  Pc 


(2.11) 


b  =  0.08664  R  Tc  /  Pc 


(2.12) 


For  hydrocarbons  and  certain  nonhydrocarbons,  the  additional 
constant  a  is  a  function  of  the  temperature  and  the  acentric 


factor : 


a  -  {1  +  S  (  1  -  /Tr  )  1 


(2.13) 


The  S  parameter  was  given  by  Graboski  and  Daubert  (17)  as 


S  =  0.48508  +  1  .5517  cj  -  0.15613  or 


(2.14) 


The  equation  of  state  was  extended  to  mixtures  by  Soave 
(67)  using  the  combining  rules  given  as: 


a  a 


ZZ  x  j  x  j  a  j  j  a  i  j 

i  j 


(2.15) 
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b  =  Z  x i b i  (2.16) 

i 

where 


aij«ij=/(3i3jai°:j)(1  -  C  j  j  )  (2.17) 

Graboski  and  Daubert  (17)  computed  the  interaction 
coefficients,  Cjj,  using  the  SRK  model,  for  a  number  of 
hydrogen  binaries.  They  restricted  the  values  of  the 
interaction  coefficients  to  C  ■,  }  <  1  and  observed  that  the 
predicted  phase  behaviour  was  insensitive  to  C  s  j  for  heavy 
hydrocarbons  at  higher  temperatures.  However,  they  found  the 
SRK  method  fitted  the  experimental  solubility  data  only  with 
unrealistically  large  positive  and  negative  values  of  the 
interaction  coefficients  for  some  of  the 
hydrogen-hydrocarbon  systems.  In  the  cases  where  Cjj  was 
greater  than  1  physically  meaningless  au  <  0  resulted.  The 
interaction  coefficients  were  also  dependent  on  temperature. 
A  similar  pattern  of  behavior  was  observed  by  El-Twaty  and 
Prausnitz  (14). 

2.3.3  Penn-State-Soave  Method 

Graboski  and  Daubert  (17)  derived  the  function  for  S  in 
equation  (2.14),  based  on  a  regression  of  vapor  pressure 
data  and  applied  the  result  to  correlate  vapor  liquid 
equilibria  data  for  hydrocarbon  mixtures  and  for  binary 
mixtures  of  hydrocarbon  with  carbon  dioxide,  hydrogen 
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sulphide,  carbon  monoxide  and  nitrogen.  The  basic  SRK 
procedure  failed  to  give  accurate  results  for  hydrogen 
K-values.  Thus,  they  proposed  the  new  expression  for  a  of 
hydrogen  as  given  in  equation  (2.18): 

a  =  1.202  exp  (-  0.30228  Tr  )  (2.18) 

This  equation  can  be  expected  to  be  accurate  for  hydrogen  at 
reduced  temperatures  greater  than  2.5.  They  proposed  that 
the  modified  a  function  (for  hydrogen)  be  used  directly  for 
hydrogen  systems  without  interaction  coefficients.  Testing 
their  modified  SRK  method,  called  the  PSS  method,  for  binary 
mixtures  of  hydrogen  in  1 -methylnaphthalene ,  bicyclohexyl, 
tetralin  and  diphenylmethane  they  found  average  deviations 
to  be  18.6%  and  10.2%  in  K-values  of  hydrogen  and  the  heavy 
solvent,  respectively. 

The  advantage  of  SRK  over  PSS  is  that  the  deviations  of 
the  calculated  results  from  experimental  values  for  the  SRK 
method  do  not  appear  to  increase  as  the  molecular  weight  of 
the  hydrocarbon  solvent  is  increased,  and  the  PSS  method 
fails  to  represent  the  K-values  of  hydrogen  in  heavy 
hydrocarbons  as  accurately  as  in  the  lighter  solvents  (29). 

In  the  case  of  hydrogen  sulphide  existing  in  the 
mixtures,  the  relationship  (2.13)  is  valid.  The  interaction 
coefficient  for  hydrogen  sulphide  and  the  hydrocarbon  is 
then  given  by  equation  (2.19): 
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Cs j (H2S)  =  0.0178  +  0.0244  A6 


(2.19) 


where 


A8  =  | 6  -  6 

HC  H2S 


(2.20) 


6  is  the  Scatchard-Hildebrand  solubility  parameter. 

2.3.4  Peng-Robinson 

The  equation  of  state  proposed  by  Peng  and  Robinson 
(41)  is  of  the  form 

RT  a 

P  =  - -  (2.21) 

v-b  v ( v+b ) -b ( v-b ) 

Equation  (2.21)  can  be  written  as 

z3  -  ( 1 -B ) z  2  +  ( A- 3B 2 -2B ) z  -  ( AB-B  2 -B  3 )  =  0  (2.22) 


where 


A  =  a  P  /  R2  T2 


B  =  b  P  /  R  T 


(2.23) 


(2.24) 


z  =  P  v  /  R  T 


(2.25) 


29 


Equation  (2.22)  can  be  solved  analytically  to  determine  the 
compressibility  factor  z.  The  equation  yields  one  or  three 
real  roots  depending  on  the  number  of  phases  in  the  system. 
In  the  two  phase  region,  the  largest  root  is  the 
compressibility  factor  of  the  vapor  while  the  smallest 
positive  root  corresponds  to  that  of  the  liquid. 

The  fugacity  coefficient  of  any  component  in  a  mixture 
is  related  to  the  volumetric  behaviour  of  the  vapor  phase. 

P 

In  0 j  =  /  (Zj-1)  d  In  P  (2.26) 

0 

Here  Zj  =  Pv  ,  /RT  and  v,  is  the  partial  molar  volume  of  the 
ith  component.  Applying  equation  (2.21)  to  (2.26)  the 
following  expression  for  the  fugacity  coefficient  is 
obtained : 


A  2Ix i a i  j  b  j 
— — ( -  -  — )  x 

2i/2B  a  b 

(2.27) 


where 


bi 

In  0  j  =  — (z- 1 )  -  In  (z-b)  - 
b 

Z+2.414B 

In  ( - ) 

Z-2.414B 


a  =  ZZ  y i y  ja i  j 

i  j 


(2.28) 
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b  =  I  y  i  b  i 


(2.29) 


i 


=  /(a  j  a  j ) (  1-6  ,  j  ) 


(2.30) 


In  equation  (2.30)  Sjj  is  an  empirically  determined  binary 
interaction  coefficient  characterising  the  binary  formed  by 
the  i  and  j  component.  It  is  equivalent  to  the  interaction 
coefficient  C|j  in  the  SRK  equation  of  state. 

The  constants  a  and  b  are  related  to  the  critical 
properties  and  the  acentric  factor  of  component  i ,  such  that 


a  j  ( T )  =  a  i  ( Tc  )  a  j  ( Tr  ,  cj  ) 


(2.31 ) 


b i  (T)  =  b  j  (Tc ) 


(2.32) 


a  i  (Tc  )  =  0.45724  R2  Tc  2  /  Pc 


(2.33) 


bj  (Tc  )  =  0.07780  R  Tc  /  Pc 


(2.34) 


a 


{  1  +  K  (  1  -  i/T,  )  } 


(2.35) 


K 


0.37464  +  1.54226  u  ~  0.26992  cj 


(2.36) 


The  fugacity  coefficient  <t>  ■,  is  defined  by  : 


<t>  j  =  f  i  /  x  j  P 


(2.37) 
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where  f-,  is  the  fugacity.  K-values  are  found  from  the 
relation 


(2.38) 


Kj  =  0s(liquid)  /  0j(vapor) 


2.3.5  Modified  Soave-Redl ich-Kwong 

El-Twaty  and  Prausnitz  (14)  found  that  when  equation 
(2.16)  was  applied  to  mixtures  of  hydrogen  and  heavy 
hydrocarbons,  the  experimental  data  could  be  fitted  only 
with  unrealistic  values  of  C  j  j .  In  some  cases,  it  was 
necessary  to  use  Cjj  >  1,  giving  meaningless  au  <0.  It  was 
shown  that  calculated  vapor-liquid  equilibria  for 
hydrogen-containing  systems  were  more  sensitive  to  small 
changes  to  b  than  to  large  changes  to  a  .  Equation  (2. 16) 
was  replaced  by  equation  (2.39): 


(2.39) 


b=IXjbi  +  xH 
i 


where  subscript  H  stands  for  hydrogen  and  the  summation  over 
i  extends  to  all  components  (including  hydrogen) ,  and  the 
summation  over  j  extends  to  all  components  except  hydrogen. 
When  EHj  =  0.,  equation  (2.39)  reduces  to  equation  (2.16). 
Using  equations  (2.9),  (2.15),  and  (2.39)  to  find  the 
fugacity  coefficient  for  component  i,  we  obtain: 
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s  1  a  v+b  b 

In  0|  =  +  Xi( -  ~  - (  In  -  -  - )}  (2.40) 

v-b  RTb2  v  v+b 


where  the  fugacity  coefficient  <j>  ■,  is  defined  by  equation 
(2.37).  0s j  designates  the  fugacity  coefficient  derived  from 
equation  (2.10).  For  hydrogen, 

XH  =  ( 1-xH)  Z  XjEHj  (2.41  ) 

j 


For  all  other  components 


X  i  =  xH  ( E H j  -  Z  x k E Hk  ) 

k 


(2.42) 


where  the  summation  over  components  k  includes  all 
components  (EHk  =  0  when  k  is  hydrogen). 

For  temperatures  above  the  critical,  the  characteristic 
constant,  a,  is  calculated  from  equation  (2.43): 

a  =  exp  {2  S(1  -  i/Tr  )}  (2.43) 

Otherwise,  equation  (2.13)  applies. 

2.4  Characterization  of  Mixtures 

Table  1  lists  the  properties  needed  to  calculate  the 
VLE  data.  The  PR,  SRK  and  Gray son-St reed  procedures  require 


• 

TABLE  1 


PROPERTIES  OF  COMPONENTS  REQUIRED  FOR  PR,  SRK 
AND  GRAYSON  -  STREED  METHODS 

PR  SRK  GS 


Critical  Temperature 

yes 

yes 

yes 

Critical  Pressure 

yes 

yes 

yes 

Acentric  Factor 

yes 

yes 

yes 

Solubility  Parameter 

no 

no 

yes 

Liquid  Molar  Volume 

no 

no 

yes 

Interaction  Parameter 

yes 

yes 

no 
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that  each  component  of  the  mixture  be  characterised  by  its 
critical  pressure,  critical  temperature,  and  acentric 
factor.  In  addition,  the  solubility  parameter  and  the  liquid 
molar  volume  of  each  component  are  used  in  the  'regular' 
solution  model  of  Chao-Seader  and  Grayson-St reed .  The  values 
of  these  parameters  for  well-defined  hydrocarbons  are 
compiled  in  Reid  et  al.  (50). 

Heavier  hydrocarbons  such  as  bitumens  and  gas  oils  are 
not  pure  compounds  but  contain  numerous  constituents  that 
would  be  difficult,  if  not  impossible,  to  identify. 
Consequently  these  hydrocarbons  are  treated  either  as: 

(i)  single  pseudo-substances, 

(ii)  fractions  derived  from  distillation  cuts,  or 

(iii)  fractions  that  are  further  divided  into  three 
different  groups,  viz.,  paraffins,  naphthenes, 
and  aromatics  (PNA) . 

The  GPA-Peng-Robinson  Programs  (43)  have  a  routine  that 
calculates  the  critical  properties  and  the  acentric  factor 
of  heptane-plus  fractions  by  the  PNA  analysis.  This  method 
requires  the  molecular  weight  (or  the  carbon  number) ,  weight 
average  or  volume  average  boiling  point  and  the  specific 
gravity  of  each  undefined  fraction  as  inputs.  The  method  of 
Cavett  (3),  and  that  described  in  the  API  Data  Book  (1)  are 
commonly  used  to  characterise  petroleum  fractions.  In  these 
methods  the  boiling  point  parameters  and  the  specific 
gravities  are  specified.  In  addition,  molecular  weights  are 
required  in  the  experimental  determination  of  molar  liquid 
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phase  compositions.  Detailed  descriptions  of  these  three 
methods  follow. 

2.4.1  Cavett  Method 

The  pseudoc r i t ical  temperature  Tc  (in  degrees  Rankine) 
is  calculated  from  equation  (2.44): 

Tc  =  a0  +  aiT  +  a2T2  +  a3AT  +  a„T3 

+  a  5 AT2  +  a  6 A 2  T 2  (2.44) 

where  T  =  molal  average  boiling  point,  °F 
A  =  0  API 

Constants  a0,  a1;...,  a6  are  listed  in  Cavett  (3). 

The  pseudocr itical  pressure  Pc  (in  psia)  is  calculated 
using  equation  (2.45): 

log  Pc  =  b0  +  biT  +  b2T2  +  b3 AT  +  b4T3 

+  b  5 AT 2  +  b6A2T  +  b7 A  2T2  (2.45) 

where  T  =  mean  average  boiling  point,  °F 

Constants  b0 ,  b1f...,  b7  are  listed  in  Cavett  (3). 

The  acentric  factor  cj  is  derived  from  the  following 

relationship: 


-  1.0 


CJ 


log  10  Pr 


(2.46) 
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where  Pr  is  reduced  vapor  pressure  at  reduced  temperature  of 
0.7. 

In  the  absence  of  vapor  pressure  data,  equation  (2.47) 
proposed  by  Edmister  (13)  is  used: 

3  In  Pc 

u  =  -  -  1  (2.47) 

7(TC/T  -1  ) 

where  Tc  and  Pc  are  calculated  from  equations  (2.44)  and 
(2.45),  respectively.  T  is  the  mean  average  boiling  point. 

2.4.2  API  Data  Book  Method 

The  boiling  point  data  of  each  fraction  is  converted  to 
ASTM  D  86  temperatures  using  Figure  3A1.1  of  Technical  Data 
Book-Petroleum  Refining  (1).  The  molal  average  boiling  point 
and  mean  average  boiling  point  are  obtained  from  Figure 
2B1.1.  The  pseudocr itical  temperature  and  pressure  are 
obtained  from  Figure  4D4 . 1 .  Figure  2B2.1,  redrawn  from  the 
Winn  nomograph,  is  used  to  calculate  the  molecular  weight, 
given  specific  gravity  and  mean  average  boiling  point.  The 
acentric  factor  is  obtained  from  Figure  2B2.3.  The 
solubility  parameter  is  derived  from  chart  developed  by 
Kesler  and  included  in  Chap.  8  of  the  Data  Book. 

The  Winn  nomograph  was  represented  analytically  by  Sim 


and  Dauber t  ( 60 ) . 
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2.4.3  PNA  Analysis 

In  this  method  it  is  assumed  that  the  fractions  can  be 
divided  into  three  groups,  namely,  n~ alkanes, 

r?-alkyl-cyclohexanes ,  and  n-alkyl-benzenes ,  based  on  their 
molecular  structures. 

The  normal  boiling  point  temperatures  have  been  fitted 
to  an  equation  of  the  form: 

In  Tb  =  Z  Cj  (n-6) j-'  (2.48) 

I 

for  each  of  the  hydrocarbon  groups.  Tb  is  the  normal  boiling 
point  temperature  (K)  and  n  is  the  number  of  carbon  atoms  in 
the  molecule.  Constants  C|  are  listed  in  Table  1  of  Peng  and 
Robinson  ( 42 ) . 

To  generate  the  critical  pressure  values  for  the 
fractions  a  group  contribution  correlation  of  Lydersen  (32) 
is  used: 


1  4 . 026n+y 

Pc  =  - 

(  0 . 227n  +  /S )  2 


(2.49) 


where  Pc  is  in  atmospheres.  The  coefficients  /3  and  y  are 
listed  in  Table  1  of  (42). 

Assuming  a  linear  function  between  acentric  factor  and 
carbon  number  the  following  equation  is  derived. 


Cl) 


f  n  +  g 


(2.50) 


. 
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where  f  and  g,  coefficients  for  all  three  groups,  are  listed 
in  Table  1  of  ( 42 ) . 

The  critical  temperature  Tc  is  calculated  from 
Edmister's  equation  (13),  which  rearranged  gives: 


3  In  Pc 

Tc  =  S  Tb  {1  +  - }  (2.51) 

7( 1+w) 


where  S  is  a  correction  factor  given  by  equation  (2.52): 


S  =  0.996704  +  0.00043155  n  (2.52) 


The  overall  parameters  for  each  cut  can  be  calculated 
according  to  the  following  equations: 


Tc  =  L  x,  Tci 
i 

Pc  =  Z  X  j  Pc  i 

i 

-  (  1  +Cl>  j  ) 

L  x iPc i 10 

i 

cj  =  -  log  { - - -  "  1 } 

Pc 


(2.53) 

(2.54) 


(2.55) 


where  subscript  i  stands  for  paraffin,  naphthenic,  and 
aromatic . 

In  the  GPA-Peng-Robinson  Programs  (43),  the  following 
options  are  permitted  in  order  to  characterise  the  undefined 
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components : 

1.  average  carbon  number  (or  molecular  weight),  specific 
gravity  and  boiling  point  for  the  entire  fraction  are 
spec i f ied , 

2.  carbon  number  (or  molecular  weight)  and  the  PNA  analysis 
are  specified, 

3.  carbon  number,  weight  average  boiling  point,  and 
molecular  weight  are  specified, 

4.  carbon  number  (or  molecular  weight),  volume  average 
boiling  point,  and  specific  gravity  are  provided,  or 

5.  similar  to  (ii)  except  that  the  PNA  analysis  is  given  in 
terms  of  volume  fractions  instead  of  mole  fractions. 

Other  methods  exist  that  characterise  particular  groups 
of  hydrocarbons  and  are  restricted  to  fractions  whose 
molecular  weights  and  boiling  points  are  far  below  those  oi 
bitumens.  Recently  Wilson  et  al.  (71)  proposed  correlations 
predicting  the  critical  properties  of  aromatic  fractions. 

It  has  become  evident  from  recent  studies  that  the 
correlation  of  hydrogen-hydrocarbon  systems  is  not  handlec 
well  by  the  SRK  equation  of  state,  particularly  at 
conditions  of  high  temperatures  and  pressures.  The  moditied 
SRK  and  the  PR  equations  were  used  to  correlate  the 
experimental  hydrogen-heavy  hydrocarbon  VLE  data.  The 

ed  with  the  Grayson-St reed  predictions. 


results  were  compar 


‘  • 
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3.  Solubility  Experiments 


3.1  Reagents  and  Materials 

The  samples  of  bitumens  and  hydrocracked  fractions  used 
in  this  study  were  obtained  from  the  Department  of  Energy, 
Mines  and  Resources,  and  were  received  with  the  following 
notations : 

1.  1-CL-77  Cold  Lake  Bitumen 

2.  4-AB-77  Athabasca  Bitumen 

3.  3-LL-77  Lloydminster  Residue 

4.  2-AB-77  650°F+  Topped  Athabasca  Bitumen 

5.  36-1  800°F+  Topped  Athabasca  Bitumen 

6.  37-1  650°F+  Topped  Cold  Lake  Bitumen 

7.  6-CL-77  800  °F+  Topped  Cold  Lake  Bitumen 

8.  3-LL-78  650  0  F+  Lloydminster  Residue 

9.  1-LL-78  7  60  0  F+  Lloydminster  Residue 

10.  4-AM-78  Amoco  Feed 

11.  91-1-1  60%  Conversion  Heavy  End 

12.  78-T-22  74.7%  Conversion  Heavy  End 

13.  78-T-23  80.1%  Conversion  Heavy  End 

14.  78-T-24  84.3%  Conversion  Heavy  End 

15.  EMR  Gas  Oil 

Samples  1  to  10  are  the  feedstocks  of  the  CANMET 
Process.  The  Heavy  Ends  (samples  1  1-14)  and  the  EMR  Gas  Oil 
are  the  products.  The  properties  of  the  solvents  1-10, 
supplied  by  EMR,  are  listed  in  Tables  2  to  4  at  the  end  of 
this  section,  along  with  those  of  a  conventional  crude  oil 
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(47).  Molecular  weights  of  solvents  1,  2,  3,  5  and  14, 
determined  by  the  Hydrocarbon  Research  Centre  using  the 
osmometric  method,  are  given  in  Table  8. 

The  boiling  point  distributions  of  the  solvents  4,  8, 
and  10  (adjusted  to  ASTM  D  2892  atmospheric  equivalent)  are 
listed  in  Tables  5-7.  The  original  Hempel  distillation  data 
for  solvents  4,  8  and  10,  supplied  by  EMR ,  are  listed  in 
Tables  48-50  of  Appendix  III.  These  data  were  not  suitable 
for  characterization  of  the  bitumen  feedstocks  due  to  their 
very  high  residue  contents. 

Mesitylene  was  obtained  from  Eastman  Organic  Chemicals 
and  had  a  boiling  range  of  160-163°C.  Hydrogen  was  obtained 
from  Linde  in  6000  psi  cylinders  with  reported  purity  of 
99.995%,  and  hydrogen  sulphide  supplied  by  Matheson  Gas 
Products  was  C.P  grade  with  minimum  purity  of  99.6%. 
Sub-bituminous  and  Lignite  coal  samples  (-200  mesh)  were 
supplied  by  the  Department  of  Energy,  Mines  and  Resources, 
Ottawa . 

3.2  Experimental  Procedure 

The  equilibrium  cell  used  was  a  2  litre,  316  stainless 
steel  Magnedrive  autoclave  (Model  AF  2005)  from  Autoclave 
Engineers  Inc.  Equilibrium  pressures  were  monitored  with 
0-340  bar  and  0-1000  psi  316  stainless  steel  Heise  Bourdon 
tube  gauges  which  had  been  calibrated  against  a  dead  weight 
gauge.  The  temperature  of  the  cell  was  measured  with  a  type 
J  iron-Constantan  thermocouple  inserted  into  a  thermowell 
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TABLE  2 


Properties  of  Cold  Lake  Bitumen  and  Conventional  Oil 


Cold  Lake  Bitumen  Conv. 

(1-CL-77)  (37-1)  (6-CL-77)  Oil 


Sp.  Gravity 

1 5/1 5  °C 

0.999 

1  .020 

1  .0333 

o 

• 

03 

1 

O 

• 

KD 

Sulphur 

wt  % 

4.30 

4.80 

5.54 

O 

• 

CM 

1 

• 

O 

Ash 

wt  % 

0.03 

0.076 

0.073 

nil 

OCR 

wt  % 

14.6 

14.4 

18.2 

1-2 

Vanadium 

ppm  (wt) 

130 

133 

153 

2-10 

Nickel 

ppm  (wt) 

92 

- 

- 

total  metals 

Carbon 

wt  % 

82.98 

83.7 

83.99 

86 

Hydrogen 

wt  % 

10.05 

10.34 

10.25 

13.5 

Nitrogen 

wt  % 

0.34 

0.46 

0.48 

0.2 

Viscosity  at  99 °C  cSt 

94 

309 

2009 

3-7 

Pitch 

wt  % 

48 

- 

- 

- 

H/C 

0.121 

0.124 

0.122 

0.157 

. 

• 

■ 

- 
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TABLE  3 


Properties  of  Athabasca  Bitumen 


(4-AB-77) 

(2-AB-77) 

(36-1 ) 

Sp.  Gravity 

1 5/1 5  °C 

1.009 

1  .023 

1  .042 

Sulphur 

wt  % 

4.48 

4.97 

5.49 

Ash 

wt  % 

0.59 

0.66 

0.73 

CCR 

wt  % 

13.3 

15.0 

19.6 

Vanadium 

ppm  (wt) 

213 

184 

193 

Nickel 

ppm  (wt) 

67 

- 

- 

Carbon 

wt  % 

83.36 

83.32 

82.49 

Hydrogen 

wt  % 

10.52 

10.27 

9.85 

Nitrogen 

wt  % 

0.43 

0.44 

- 

Viscosity  at  99 0 C  cSt 

210 

413 

2237 

Pitch 

wt  % 

51  .5 

- 

- 

H/C 

0.126 

0.123 

0.119 
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TABLE  4 


Properties  of  Lloydminster  Residue  and  Amoco  Feed 


Lloydminster  Residue 

Amoco 

• 

(3-LL-77) 

(3-LL-78) 

( 1-LL-78) 

Feed 

Sp.  Gravity 

1 5/1 5  °C 

1  .033 

1.000 

1  .032 

0.95 

Sulphur 

wt  % 

5.02 

4.34 

5.30 

3.42 

Ash 

wt  % 

0.04 

0.035 

0.067 

0.  1 

Vanadium 

ppm  (wt) 

124 

146 

210 

124 

Nickel 

ppm  (wt) 

— 

66 

92 

37 

Carbon 

wt  % 

83.71 

86.16 

82.57 

82.25 

Hydrogen 

wt  % 

9.91 

10.57 

10.11 

11.37 

Nitrogen 

wt  % 

0.34 

0.31 

0.63 

0.29 

Viscosity  a 

t  99 °C  cSt 

3691 

188 

— 

12 

H/C 

0.118 

0.123 

0.122 

0.138 

- 
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TABLE  5 

Boiling  Point  Distribution  for  Topped 
Athabasca  Bitumen  (2-AB-77) 

Temperature  °C 


IBP- 

175- 

200- 

225- 

250- 

275- 

307- 

336- 

364- 

392- 

419+ 

175 

200 

225 

250 

275 

307 

336 

364 

392 

419 

%  Vol 

0.3 

0.4 

0.6 

1 .2 

6.2 

1.7 

1  .  1 

3.9 

7 . 5 

12.2 

64.6 

Sp  Gr 

.824 

.824 

.824 

.838 

.868 

.  88  1 

.902 

.  927 

.940 

.957 

1  .058 

S  % 

3.08 

3.08 

3.08 

2 . 83 

2.84 

2.93 

3. 10 

2 . 98 

3.03 

3.23 

5.09 

N  % 

.006 

.006 

.006 

.008 

.013 

.018 

.024 

.062 

.089 

.  132 

.667 

IBP 

1  52  0  C 

TABLE  6 

Boiling 

Point 

Distribution 

for 

L 1 oydm 

i  nster 

Res i due 

(3-LL 

-78) 

T  emperature 

°C 

IBP- 

225- 

250- 

275- 

307- 

336- 

364- 

392- 

419  + 

225 

250 

275 

307 

336 

364 

392 

419 

%  Vol 

0.3 

1 . 2 

2.0 

1 .2 

3.9 

6 . 3 

5.7 

9 . 2 

69 . 3 

IBP  2 1 8  0  C 


TABLE  7 


Boiling  Point  Distribution  for  Amoco  Feed  (4-AM-78) 

Temperature  °C 


IBP- 

200- 

225- 

250- 

275- 

307- 

336- 

364- 

392- 

419+ 

200 

225 

250 

275 

307 

336 

364 

392 

4  19 

%  Vol 

3 . 3 

8.0 

9.9 

10.4 

7 . 7 

7.2 

3.9 

3 . 7 

6 . 5 

38 . 1 

Sp  Gr 

.800 

.820 

.  843 

.  859 

.  874 

.  889 

.912 

.935 

.964 

1  .065 

S  % 

<0.3 

<0.3 

<0.3 

0.41 

0.63 

0.92 

1 .57 

2.42 

3.42 

5 . 98 

N  % 

trace 

trace 

.002 

.004 

.010 

.017 

.037 

.073 

.  130 

.  710 

IBP  1 76 0  C 


€ 


. 


. 

s 

TABLE  8 


Molecular  Weights  of  Selected  Solvents 


1-CL-77  Cold  Lake  Bitumen  515 
4-AB-77  Athabasca  Bitumen  522 
3-LL-77  Lloydminster  Residue  864 
36-1  800°F+  Topped  Athabasca  Bitumen  732 
78-T-24  84.3%  Conversion  Heavy  End  332 
EMR  Gas  Oil  250 1 


1  by  API  Method 
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extending  into  the  liquid  phase.  The  thermocouple  was 
calibrated  against  a  platinum  resistance  thermometer  for  use 
in  the  range  0°  to  402°C.  Temperatures  were  controlled  to 
within  ±0.5°C.  The  details  of  the  apparatus  are  shown  in 
Figure  2. 

For  each  run,  approximately  500  cc  of  solvent  was 
charged  into  the  autoclave  and  the  system  flushed  with 
hydrogen.  The  reactor  was  brought  to  the  desired  temperature 
and  the  hydrogen  pressure  was  set  approximately  by  adjusting 
the  cylinder  delivery  pressure.  Stirring  was  commenced  and 
maintained  during  the  entire  process  at  600  rpm,  except 
during  sample  withdrawal.  Pressure  declined  as  a  result  of 
hydrogen  dissolution.  The  cell  was  left  overnight  to  ensure 
equilibrium  was  achieved. 

Samples  of  the  vapor  were  withdrawn  from  the  top  of 
autoclave  and  analysed  using  a  Hewlett-Packard  Model  5710A 
gas  chromatograph  with  a  1/4  inch  O.D  2  metre  long  column 
containing  Porapak  QS ,  50-80  mesh  packing.  The  column  was 
operated  isothermally  at  90°C  with  a  high  purity  nitrogen 
carrier  gas  flowing  at  25  ml/min.  A  thermal  conductivity 
detector  was  used.  Typical  vapor  analysis  showed  hydrogen 
concentration  in  excess  of  99%  by  volume. 

The  stirrer  was  stopped  before  liquid  samples  were 
withdrawn  from  the  cell  into  a  weighed  sampling  flask.  The 
gases  evolved  from  the  liquid  were  measured  volumet r ical ly 
by  displacement  of  mercury  from  a  calibrated  burette.  The 
amount  of  gas  remaining  in  the  bitumen  at  atmospheric 


Legend 
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Figure  2  Schematic  Diagram  of  the  Apparatus 
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temperature  and  pressure  was,  at  most,  of  the  order  of  1%  of 
the  total  gas  evolved.  After  volume  measurement  the  gases 
were  passed  through  a  gas  chromatograph  for  analysis.  From 
this  analysis,  knowing  the  barometric  pressure,  ambient 
temperature,  and  total  gas  volume,  the  weight  of  hydrogen 
dissolved  in  the  solvent  was  computed  using  the  ideal  gas 
law.  By  weighing  the  loaded  sampling  flask,  the  weight  of 
the  solvent  was  determined  by  difference,  and  the  solubility 
calculated  as  grams  hydrogen  per  gram  solvent.  From  the 
molecular  weights  of  the  solvent  and  the  solute  gas,  the 
solubility  was  converted  to  mole  fraction. 

A  similar  procedure  was  used  to  determine  the 
solubility  of  hydrogen  sulphide  and  carbon  dioxide  in  these 
solvents.  In  the  measurement  of  hydrogen  solubility  in 
slurries  the  stirrer  was  only  slowed  down  during  sampling  to 
prevent  coal  particles  from  settling. 


3 


4.  Results  and  Discussion 

To  confirm  the  sampling  and  analytical  procedures,  the 
solubility  of  hydrogen  in  mesitylene  was  determined  at 
200  °C .  The  results  are  presented  in  Table  12  of  Appendix  I 
and  are  compared  with  data  available  in  the  literature 
(12,44)  at  204°C  in  Figure  3.  The  data  of  Eakin  and  DeVaney 
(12)  are  for  a  ternary  mixture  containing  hydrogen  sulphide. 
The  present  results  are  in  agreement  with  the  literature 
values  within  a  few  percent. 

4.1  Solubility  of  Hydrogen  in  Bitumens 

Experimental  data  for  the  solubility  of  hydrogen  in 
bitumens,  hydrocracked  products  and  EMR  Gas  Oil,  reported  as 
g  hydrogen/g  bitumen,  are  listed  in  Tables  13  to  27  of 
Appendix  I.  The  results  are  plotted  on  Figures  4  to  15.  Most 
of  the  data  were  obtained  at  200°  and  300°C;  however  a  few 
points  were  measured  at  100°,  250°  and  370°C.  For  a 

particular  sample,  at  a  given  pressure,  the  solubility  of 
hydrogen  increases  with  increasing  temperature,  an  effect 
typical  of  hydrogen  systems. 

Comparisons  of  the  solubility  plots  leads  to  the 
following  observations: 

1.  The  solubility  of  hydrogen  in  Athabasca  and  Cold  Lake 
bitumens  is  essentially  the  same. 

2.  Results  for  the  650°F+  and  800°F+  topped  bitumen 
products,  with  the  exception  of  the  800°f+  Topped  Cold 
Lake  Bitumen  (6-CL-77)  at  200°C,  support  a  conclusion 
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that  hydrogen  solubility  decreases  with  increasing 
average  molecular  weight  of  the  bitumen  samples.  The 
solubility  in  6-CL-77  at  200°C  is  higher  than  expected. 
The  solubility  of  hydrogen  is  essentially  the  same  for 
the  two  650°F+  products  and  for  the  two  800°F+  products 
except  for  the  Cold  Lake  800°F+  topped  product  at  200°C. 
The  solubility  in  this  topped  Cold  Lake  product  is 
approximately  15%  higher  than  that  in  the  topped 
Athabasca  bitumen  at  200°C. 

3.  The  solubility  of  hydrogen  in  Lloydminster  heavy  oil  is 
approximately  15%  less  than  in  Athabasca  bitumen. 
However  the  solubility  in  the  650°F+  and  760°F+ 
Lloydminster  residues  is  about  the  same  as  for  the 
Athabasca  Bitumen. 

4.  The  solubility  of  hydrogen  in  Amoco  feed  is 
approximately  25%  higher  than  in  Athabasca  Bitumen. 

5.  The  solubility  of  hydrogen  in  EMR  gas  oil  is 
approximately  10%  higher  than  in  Amoco  Feed. 

6.  The  solubility  of  hydrogen  in  60%  Conversion  Heavy  Ends 
is  higher  than  in  the  74.7%  to  84.3%  Conversion  Heavy 
Ends  at  hydrogen  partial  pressures  greater  than 
approximately  8.0  MPa. 

Significant  cracking  occurred  for  all  samples 
investigated  at  370°C;  thus  comparisons  between  samples  of 
data  obtained  at  this  temperature  are  not  valid. 

The  linear  relationship  between  pressure  and  solubility 
enables  one  to  estimate  an  apparent  Henry  s  constant,  Hr 


' 
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FIG  3  H2  SOLUBILITY  IN  MESITYLENE  AT  204  C 
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FIG  4  SOLUBILITY  OF  H2  IN  flTHRBflSCR  BITUMEN 
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FIG  5  SOLUBILITY  OF  H2  IN  COLD  LAKE  BITUMEN 
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FIG  6  SOLUBILITY  OF  H2  IN  LLOYDMINSTER  RESIDUE 
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FIG  7  SOLUBILITY  OF  H2  IN  650  F+  TOPPED 
ATHABASCA  BITUMEN 
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FIG  8  SOLUBILITY  OF  H2  IN  800  F+  TOPPED 

ATHABASCA  BITUMEN 
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FIG  9  SOLUBILITY  OF  H2  IN  650  F+  TOPPED 
COLD  LAKE  BITUMEN 
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FIG  10  SOLUBILITY  OF  H2  IN  800  F+  TOPPED 

COLD  LAKE  BITUMEN 
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SOLUBILITY  OF  H2  IN  LL0YDMINSTER  RESIDUE 
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FIG  12  SOLUBILITY  OF  H2  IN  AMOCO  FEED 
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FIG  13  SOLUBILITY  OF  H2  IN  607.  CONVERSION 

HEAVY  END 
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FIG  14  SOLUBILITY  OF  H2  IN  74.7%  CONVERSION 

HEAVY  END 
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FIG  15  SOLUBILITY  OF  H2  IN  84.37.  CONVERSION 

HEAVY  END 
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equal  to  the  slope  of  the  pressure-solubility  line.  The 
lower  the  value  of  H,  the  higher  is  the  solubility,  and  vice 
versa.  From  these  data,  the  effect  of  solvent  molecular 
weights  on  the  solubility  can  be  studied  by  plotting  H 
against  solvent  molecular  weight.  Table  9  lists  the  solvent 
molecular  weights  and  the  Henry’s  constants,  H,  at  200°  and 
300°C.  The  plot  of  H  versus  molecular  weight  is  shown  on 
Figure  16.  There  is  sharp  increase  in  H  up  to  molecular 
weight  of  350  and  then  a  slow  increase  above  350.  Also,  H 
increases  with  decrease  in  temperature,  as  expected. 

4.2  Solubility  of  Hydrogen  in  Slurries 

The  results  for  the  solubility  of  hydrogen  at  200°  and 
300°C  in  a  10  wt.%  sub-bituminous  coal-90  wt.%  Cold  Lake 
bitumen  slurry,  reported  as  g  hydrogen/g  bitumen  are  given 
in  Table  28.  The  results  are  plotted  on  Figure  17.  The 
solubility  is  essentially  the  same  as  in  the  Cold  Lake 
Bitumen.  Results  for  40  wt.%  slurries  of  Cold  Lake  Bitumen 
and  sub-bituminous  and  lignite  coals  respectively  are  given 
in  Tables  29  and  30  and  Figures  18  and  19  for  temperatures 
of  200°,  250°  and  300°C.  The  results  are  interesting  and 
show  a  decrease  in  solubility  with  an  increase  in 
temperature  up  to  a  pressure  of  approximately  16.5  MPa  and 
then  an  increasing  solubility  with  increasing  temperature  at 


. 
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TABLE  9 


Effect  of  Molecular  Weight  on  Apparent  Henry’s  Constant 


Mol  Wt 

Apparent  ] 

Henry's  Constant 

( 

MPa  (kg 

bitumen/g  H2)  ) 

200  °C 

300  0  C 

250 

10.78 

7.59 

332 

15.31 

10.85 

515 

16.22 

11.50 

522 

16.50 

11.73 

864 

17.47 

12.53 

( 

MPa  (kg 

bitumen/g  H2S ) ) 

200  0 C 

300  0 C 

332 

0.106 

0.297 

522 

0.134 

0.218 

( 

MPa  (kg 

bitumen/g  C02 ) ) 

50  °C 

1  50  0  C 

250 

0.0638 

0.1312 

50  °C 

300  0 C 

522 

0.0865 

0.2052 
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FIG  16  EFFECT  OF  SOLVENT  MOLECULAR  WEIGHT 
ON  APPARENT  HENRY’S  CONSTANT 
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FIG  17  HYDROGEN  SOLUBILITY  IN  10  WT  / 
SUB-BITUMINOUS  COAL/BITUMEN  SLURRY 
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FIG  18  HYDROGEN  SOLUBILITY  IN  40  WT  / 
SUB-BITUMINOUS  COflL/BITUMEN  SLURRY 
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FIG  19  HYDROGEN  SOLUBILITY  IN  40  WT  7. 
LIGNITE  CORL/BITUMEN  SLURRY 
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FIG  20  HYDROGEN  SOLUBILITY  IN  25  WT  7. 
SUB-BITUMINOUS  COflL/BITUMEN  SLURRY 
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higher  pressures.  Measurements  were  made  in  a  random  order 
and  with  several  different  initial  charges  to  confirm  this 
result.  The  type  of  coal  does  not  appear  to  have  a 
significant  effect  on  the  hydrogen  solubility.  The 
solubility  is  higher  in  the  slurry  than  in  Cold  Lake  Bitumen 
at  200°C  but  lower  at  300°C  at  hydrogen  partial  pressures  up 
to  20.0  MPa.  At  higher  partial  pressures,  at  300°C,  the 
solubility  appears  to  be  the  same  as  in  Cold  Lake  Bitumen. 

Data  for  a  25  wt . %  slurry  are  given  in  Table  31  and  are 
plotted  on  Figure  20.  The  solubility  of  hydrogen  for  this 
slurry  increases  with  increasing  temperature  over  the  range 
of  hydrogen  partial  pressures  investigated  and  is  higher 
than  in  Cold  Lake  bitumen. 

4.3  Solubility  of  H2S  and  C02  in  Bitumens 

Equilibrium  solubility  data  for  hydrogen  sulphide  in 
Athabasca  Bitumen  and  Heavy  End  (78-T-24)  are  given  in 
Tables  32  and  33  and  are  plotted  on  Figures  21  and  22.  The 
solubility  decreases  with  rising  temperatures  and  is  the 
same  for  both  samples  at  200°C.  At  300°C,  the  solubility  of 
hydrogen  sulphide  is  higher  in  the  Athabasca  Bitumen. 
Solubility  data  for  the  EMR  Gas  Oil  are  given  in  Table  34 
for  100°,  200°  and  300°C.  The  data  for  100°  and  300°C  are 
plotted  on  Figure  23.  The  solubility  is  higher  in  the  EMR 
Gas  Oil  than  in  the  heavier  Athabasca  Bitumen. 

The  apparent  Henry's  constants,  H,  for  hydrogen 
sulphide  are  listed  in  Table  9.  They  are  estimated  at  zero 
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solubilities.  The  values  increase  with  increasing 
temperature . 

Carbon  dioxide  solubility  data  in  Athabasca  Bitumen  and 
EMR  Gas  Oil  are  listed  in  Tables  35  and  36,  respectively. 
Temperatures  range  from  50°  to  200°C,  with  pressures  up  to 
6.4  MPa.  The  data  are  plotted  on  Figures  24  and  25.  The 
solubility  decreases  with  rising  temperatures,  as  expected. 
It  is  higher  in  the  EMR  Gas  Oil  compared  to  that  in  the 
higher  molecular  weight  Athabasca  Bitumen. 

The  apparent  Henry's  constants,  H,  are  listed  in  Table 
9  for  carbon  dioxide  in  EMR  Gas  Oil  at  50°  and  150°C  and  in 
Athabasca  Bitumen  at  50°  and  300°C.  The  values  increase  with 
increasing  temperature  and  increasing  solvent  molecular 
weight . 
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FIG  21  SOLUBILITY  OF  H2S  IN  ATHABASCA 

BITUMEN 


H2S  PRESSURE.  MPfi 
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FIG  22  SOLUBILITY  OF  H2S  IN  84.3% 
CONVERSION  HEAVY  END 


H2S  PRESSURE,  MPR 
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FIG  23  SOLUBILITY  OF  H2S  IN  EUR  GRS  OIL 
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FIG  24  SOLUBILITY  OF  C02  IN  ATHABASCA 

BITUMEN 
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FIG  25  SOLUBILITY  OF  C02  IN  EMR  GAS  OIL 


5.  Correlation  and  Fraction  Characterization 

One  of  the  useful  applications  of  the  experimental  data 
is  to  test  the  validity  of  existing  predictive  methods,  such 
as  the  methods  of  Chao-Seader,  Grayson-St reed , 
Soave-Redl ich-Kwong ,  modified  Soave-Redl ich-Kwong  and 
Peng-Robinson .  In  practice,  to  do  VLE  calculations  by  any  of 
these  equations,  the  following  procedures  are  implemented. 
An  initial  set  of  K-values  are  assumed  and  a  flash 
calculation  is  performed  to  find  the  liquid  and  vapor 
compositions  of  each  component.  Using  the  calculated  phase 
compositions,  the  equation  of  state  can  be  solved  for  the 
fugacities  of  the  liquid  phase  and  the  vapor  phase.  If  the 
fugacities  are  not  equal  for  all  the  components  in  both  the 
phases,  updated  K-values  are  used  and  the  flash  calculation 
repeated  until  the  fugacities  match. 


5.1  Results  and  Discussion 


The  presence 

of  a 

light  gas,  like 

hydrogen , 

in 

hydrocarbon 

systems  is 

not  handled  very 

well  by 

the 

Chao-Seader 

and 

the 

SRK  correlations 

at  eleva 

ted 

temperatures  and  pressures.  The  high  temperatures  required 
in  the  hydrotreating  processes  exceed  the  upper  limit  value 
of  260  0 C  set  for  hydrogen  in  the  Chao-Seader  method.  To 
predict  the  VLE  behavior  of  hydrogen-hydrocarbon  systems  by 
the  SRK  equation  requires  large  interaction  coefficients, 
Cjj.  This,  however,  leads  to  physically  meaningless  values 
in  the  a  term. 
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It  is  not  immediately  evident  that  an  equation  of  state 
can  successfully  predict  VLE  properties  of  hydrogen-heavy 
hydrocarbon  systems.  The  PR,  modified  SRK  and  the 
Grayson-Streed  methods  were  used  in  this  work  to  correlate 
experimental  VLE  data  over  wide  ranges  of  temperature, 
pressure,  and  molecular  weight  of  solvents.  Toward  that  end, 
the  experimental  vapor  and  liquid  compositions  reported  by 
Sebastian  et  al.  (53,55,56),  Simnick  et  al.  (61,62,64),  Yao 
et  al.  (72)  and  Nichols  et  al.  (39)  for  hydrogen  in  heavy 
naphthenic  and  aromatic  solvents,  along  with  data  obtained 
in  this  work  for  the  EMR  Gas  Oil  and  the  bitumens  were 
compared  with  the  values  calculated  by  the  three  methods. 
Typical  VLE  calculation  results  are  given  in  Tables  37-42  of 
Appendix  1 1 . 

The  experimental  K-values  of  hydrogen  in  bitumens  and 
EMR  Gas  Oil  were  approximated  by  1/x  (x  is  the  hydrogen 
liquid  phase  composition),  since  the  vapor  phase  was  almost 
all  hydrogen.  The  calculated  K-values  were  derived  from 
equation  (2.5)  for  Grayson-Streed  and  equation  (2.38)  for 
the  PR  and  modified  SRK  correlations.  The  critical 
properties  and  the  acentric  factors  of  the  solvents  used  in 
the  calculations  were  those  reported  in  Reid  et  al.  (50). 

5.1.1  Modified  Soave-Redlich-Kwong  Method 

The  modified  SRK  method  was  used  to  correlate  the  VLE 
behavior  of  hydrogen  in  bicyclohexyl  (56),  tetralin  (61), 
quinoline  (53),  diphenylmethane  (62),  m-cresol  (64)  and 


.  I 


81 


Athabasca  Bitumen  at  experimental  temperatures  and 
pressures.  In  all  the  cases  tested,  the  characteristic 
constant,  a,  for  hydrogen  was  calculated  using  equation 
(2.43)  and  those  for  the  solvents  were  obtained  using 
equation  (2.13).  The  R-values  of  hydrogen  were  calculated 
using  equation  (2.38)  and  the  liquid  and  vapor  phase 
fugacity  coefficients  were  obtained  from  equation  (2.40). 
The  interaction  coefficient,  Ciir  was  set  at  zero. 

By  adjusting  the  binary  parameter  EHj  for  each 
isotherm,  very  good  agreement  between  the  calculated  and 
experimental  hydrogen  liquid  phase  compositions  was 
obtained.  To  determine  the  optimum  E  Hj  for  each  isotherm  the 
average  absolute  deviation  (AAD)  of  the  calculated  hydrogen 
liquid  phase  compositions  from  experimental  values  was 
computed.  The  optimum  EHj  for  an  isotherm  was  one  giving  the 
minimum  AAD.  These  are  given  in  Table  10,  and  the  average 
absolute  deviations  are  given  in  Table  11. 

Figure  26  shows  the  results  for  the  K-values  of 
hydrogen  for  hydrogen-bicyclohexyl  system  at  189°  and  349°C. 
The  E  H  j  values  used  were  7.4  and  21.4,  respectively,  giving 
AAD  of  2.2  %  and  3.3  %.  These  low  errors  show  the 
hydrogen-bicyclohexyl  system  is  well  represented  by  the 
modified  SRK  equation  of  state  based  on  equation  (2.39)  with 
C  j  j  of  zero.  The  VLE  data  are  listed  in  Table  37. 

Table  10  shows  that  the  hydrogen-heavy  hydrocarbon 
systems  are  well  handled  by  the  modified  SRK  equation. 
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10°  101  102 

PRESSURE  .  MPfl 


FIG  26  HYDROGEN  K-VALUE  VS  PRESSURE 
FOR  BICYCLOHEXYL  (SEBASTIAN  ET  AL( 56 )  ) 
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TABLE  10 


INTERACTION  PARAMETERS  REQUIRED  IN  PR, 
MODIFIED  PR  AND  MODIFIED  SRK  METHODS 


SOLVENT 


BI CYCLOHEXYL 


TETRALIN 


QUINOLINE 


r  (°c) 

*i  j 

E  H  j (Mod  PR) 

E  H j (Mod 

189.0 

-1.4 

3.5 

7.4 

268.7 

-1.4 

5.7 

348.6 

-1.4 

11.0 

21.4 

428.5 

-1.4 

-  33.0 

189.6 

-0.6 

4.6 

8.0 

268.7 

-0.6 

7.6 

13.4 

348.6 

-0.6 

13.8 

22.4 

389.  1 

-0.6 

18.4 

189.3 

-3.4 

-2.9 

268.7 

-5.0 

-3.5 

348.6 

-6.7 

-4.  1 

428.5 

-8.8 

-4.5 

189.6 

-0.8 

3.0 

5.9 

268.7 

-0.8 

3.9 

348.6 

-0.8 

6.8 

15.7 

428.5 

-0.8 

16.7 

DIPHENYLMETHANE 
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TABLE  10  Continued... 


1 -METHYLNAPHTHALENE 

189.0 

o 

• 

o 

4.0 

268.7 

o 

• 

o 

5.3 

348.6 

o 

• 

o 

9.2 

428.5 

o 

• 

o 

18.5 

THIANAPHTHENE 

187.7 

0.3 

8.3 

267.9 

0.3 

10.9 

348.0 

0.3 

17.4 

429.5 

0.3 

29.0 

/77-CRESOL 

188.9 

2.7 

4.7 

268.6 

3.8 

7.6 

348.3 

7.4 

14.5 

H-HEXANE 

4.4 

-0.6 

-1.0 

37.8 

KO 

• 

o 

1 

71.1 

• 

o 

1 

1  .5 

137.8 

-0.6 

6.0 

171.1 

<o 

• 

o 

1 

11.1 

ATHABASCA  BITUMEN 

200.0 

-4.6 

-8.6 

-4.  1 

300.0 

-4.6 

-19.9 

i 

K> 

• 

O 

COLD  LAKE  BITUMEN 

200.0 

-4.2 

-6.6 

300.0 

-4.2 

-21.1 

EMR  GAS  OIL 

200.0 

-3.2 

-2.  1 

-3.2  0.1 


300.0 


• 
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TABLE  1 1 


COMPARISON  OF  AVERAGE  ABSOLUTE  DEVIATIONS 


SOLVENT 

o 

o 

PR 

Mod  PR 

Mod  SRK 

GS 

““  “ “  “■  “■ 

—  — 

BI CYCLOHEXYL 

189 

2.75 

3.  1 

2.2 

18.2 

268.7 

8.6 

0.6 

23.9 

348.6 

10.7 

2.  1 

3.3 

9.0 

428.5 

9.2 

7.4 

12.9 

TETRALIN 

189.6 

0.57 

0.66 

1  .2 

23.9 

268.7 

1.4 

1.8 

1  .5 

32.8 

348.6 

6.2 

3.0 

4.4 

14.7 

389.  1 

6.4 

3.7 

8.7 

QUINOLINE 

189.3 

1  .6 

1.9 

268.7 

0.65 

0.85 

348.6 

1  .  1 

1  .7 

428.5 

0.36 

0.37 

DIPHENYLMETHANE 

189.6 

1.4 

0.33 

1.0 

16.9 

268.7 

8.4 

0.98 

22.9 

348.6 

7.9 

1  .3 

1.5 

20.5 

428.5 

5.6 

1  .0 

8.0 

1 -METHYLNAPHTHALENE 189.0 

1  .4 

1.7 

28.7 

268.7 

1  .  1 

1.3 

38.3 

348.6 

7.  1 

1  .  1 

4  1.2 

428.5 

12.1 

2.3 

29.9 

• 

• 
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TABLE  11  Continued... 

THIANAPHTHENE  187.7 

267.9 

348.0 

429.5 

/77-CRESOL  188.9 

268.6 
348.3 


n-HEXANE 

4.4 

2.6 

37.8 

2.7 

71.1 

3.4 

137.8 

171.1 

ATHABASCA  BITUMEN  200.0 

300.0 

COLD  LAKE  BITUMEN  200.0 

300.0 
200.0 
300.0 


0.93 

0.63 

4.4 

1.7 

14.1 

2.  1 

26.8 

6.0 

0.48 

o 

• 

00 

0.77 

1  .4 

1  .4 

2.6 

1  .4 

1.8 

2.2 

2.7 


1.7 

0.34 

13.3 

35.8 

0.41 

12.6 

0.63 

1  .4 

21.5 

37.8 

1.6 

12.7 

0.97 

0.54 

18.9 

1.5 

GAS  OIL 


■ 
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Whereas  the  small  values  of  EHj  are  physically  credible,  the 
large  values  of  C  }  required  in  the  original  SRK  are  not. 
Except  for  hydrogen  in  quinoline  and  Athabasca  Bitumen,  the 
parameters  are  positive  and  increase  with  increase  in 
temperature . 

5.1.2  Peng-Robinson  Method 

The  Peng-Robinson  method  was  used  to  predict  the  VLE 
properties  for  hydrogen  binaries  in  bicyclohexyl  (56), 
tetralin  (61),  diphenylmethane  (62),  1 -methy lnaphthalene 
(72),  thianaphthene  (55),  r?-hexane  (39),  Athabasca  Bitumen, 
Cold  Lake  Bitumen  and  EMR  Gas  Oil  at  experimental 
temperatures  and  pressures.  The  characteristic  constant,  a, 
for  hydrogen  and  the  solvents  was  calculated  from  equation 
(2.35).  The  K-values  of  hydrogen  were  calculated  using 
equation  (2.38),  and  the  liquid  and  vapor  phase  fugacity 
coefficients  were  obtained  from  equation  (2.27).  The  optimum 
binary  interaction  coefficients,  6  f  j ,  were  determined  by  the 
same  procedure  used  to  find  the  optimum  EHj  in  the  modified 
SRK  method. 

No  simple  relationship  between  interaction  coefficients 
and  temperatures  was  observed.  An  overall  interaction 
coefficient,  optimized  for  the  entire  experimental 
temperature  range,  was  determined  for  each  binary.  These  are 
given  in  Table  10  and  the  AAD  values  are  given  in  Table  11. 

When  optimum  interaction  coefficient  values  are  used  to 
calculate  VLE  data,  a  good  agreement  between  the 
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experimental  and  calculated  hydrogen  liquid  phase 
compositions  exists  at  temperatures  below  250°C.  The 
deviations  rise  sharply  above  this  temperature.  For  example, 
the  hydrogen-bicyclohexyl  system  requires  an  overall 
interaction  coefficient  of  -1.4  with  AAD  of  2.8%  at  189°C 
and  10.7%  at  349°C  (cf.  2.2%  and  3.3%  for  the  modified  SRK ) . 
The  results  for  K-value  of  hydrogen  are  shown  on  Figure  27 
at  1 89 °C  and  349°C.  The  VLE  data  are  given  in  Table  38. 
Equation  (2.38)  calculates  the  K-values  which  agree 
reasonably  well  with  experiments  at  189°C  but  overpredicts 
for  pressures  above  10  MPa  at  349°C. 

Except  for  solutions  of  hydrogen  in  bicyclohexyl, 
Athabasca  Bitumen,  Cold  Lake  Bitumen  and  EMR  Gas  Oil  the 
interaction  coefficients  lie  in  the  range  -1  <  6  s  j  <  +1.  More 
unrealistic  interaction  coefficient  values  are  needed  for 
hydrogen  binaries  in  Athabasca  Bitumen,  Cold  Lake  Bitumen 
and  EMR  Gas  Oil.  If  the  critical  properties  and  acentric 
factors  of  the  solvents  are  determined  by  the  Cavett's 
method,  the  binary  interaction  coefficients  needed  to  handle 
the  three  systems  by  the  PR  correlation  are  -4.6  for  the 
Athabasca  Bitumen,  -4.2  for  the  Cold  Lake  Bitumen  and  -3.2 
for  the  EMR  Gas  Oil.  At  200°C,  the  AAD  values  are  low,  being 
1.7%  for  Athabasca  Bitumen,  0.6%  for  Cold  Lake  Bitumen  and 
1.0%  for  EMR  Gas  Oil.  However,  at  300°C,  the  errors  are 
large  because  the  calculated  VLE  are  not  sensitive  to  large 
changes  in  the  interaction  coefficients  (a  for  hydrogen 
equals  zero  at  289°C).  The  AAD  values  for  300°C  are  36%  for 
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FIG  27  HYDROGEN  K-VALUE  VS  PRESSURE 
FOR  BICYCLOHEXYL  (SEBASTIAN  ET  AL( 56 } ) 
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Athabasca  Bitumen,  38%  for  Cold  Lake  Bitumen  and  19%  for  EMR 
Gas  Oil. 

If  the  critical  properties  and  the  acentric  factors  of 
Athabasca  Bitumen  and  Cold  Lake  Bitumen  are  determined  by 
the  PNA  analysis,  the  interaction  coefficients  are  still 
large  and  negative.  The  interaction  coefficients  are  -2.7 
for  Athabasca  Bitumen  and  -2.2  for  Cold  Lake  Bitumen.  The 
AAD  values  are  0.6%  at  200°C  and  20%  at  300°C  for  Athabasca 
Bitumen,  and  1.3%  at  200°C  and  22%  at  300°C  for  Cold  Lake 
Bitumen . 

El-Twaty  and  Prausnitz  (14)  improved  the  SRK  method  by 
modifying  the  constant  b  in  the  repulsive  term  in  the 
equation.  For  mixtures,  they  proposed  equation  (2.39): 

b  =  L  x  j  b  j  +  xH  2  x  j  E  H  j  (2.39) 

1  j 

They  also  changed  the  functional  form  of  the  dependence  of 
the  characteristic  constant,  a,  on  acentric  factor  and 
temperature  by  equation  (2.43),  whenever  T  >  Tc .  The  PR 
equivalent  of  equation  (2.43)  is  given  by  equation  (4.1): 

a  =  exp  (2  k(  1  -  /Tr  )}  (4.1) 

These  two  changes  produced  significant  improvement  in 
the  correlation  of  hydrogen-heavy  hydrocarbons  binary  data 
over  a  wide  range  of  temperatures  and  pressures;  and  at  the 
same  time  requiring  physically  reasonable  values  of  E  i . 
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This  approach  was  applied  to  the  PR  equation  of  state, 
equation  (2.21).  The  new  expression  for  the  fugacity 
coefficient  is  then  given  by  equation  (4.2): 

P  Xi  A  Z+2.414B 

In  4>-,  =  In  0i  +  — { ( z- 1  )  +  -  In  ( - )}  (4.2) 

b  2/2B  Z-2.414B 

where  <pP\  designates  the  fugacity  coefficient  found  for  the 
PR  equation  ;  it  is  given  by  equation  (2.27). 

The  modified  PR  method  was  used  to  predict  VLE  behavior 
of  hydrogen  in  bicyclohexyl  (56),  tetralin  (61),  quinoline 
(53),  diphenylme thane  (62),  1 -methynaphthalene  (72), 

thianaphthene  (  55  ),  /77-cresol  (64),  n-hexane  (39),  Athabasca 
Bitumen,  Cold  Lake  Bitumen  and  EMR  Gas  Oil.  The 
characteristic  constant,  a,  was  determined  by  equation  (4.1) 
for  hydrogen  and  equation  (2.35)  for  the  solvents.  The 
K-values  of  hydrogen  were  calculated  using  equation  (2.38) 
and  the  liquid  and  vapor  phase  fugacity  coefficients  were 
obtained  from  equation  (4.2). 

The  optimum  EHj  values  required  in  (2.39)  for  each 
binary  were  determined  by  fitting  experimental  solubility 
data,  following  the  same  criterion  to  determine  the  optimum 
E  Hj  in  the  modified  SRK  method  and  5 ;  j  in  the  PR  method. 
These  are  given  in  Table  10.  The  AAD  values  are  given  in 
Table  11.  The  binary  interaction  coefficients  were  set  equal 


to  zero. 
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The 

use  of  E  Hj 

results 

in  a  good 

correlation  of 

the 

hydrogen 

solubility  data  by  the 

modi f ied 

PR 

equation . 

The 

K-values 

of  hydrogen 

plotted 

against 

pressure  for 

the 

hydrogen- 

•bicyclohexyl 

data  of 

Sebastian 

et 

al.  (56) 

are 

shown  on  Figures  26  and  27  for  comparison  with  the  modified 
SRK  and  the  original  PR  methods  (the  VLE  data  are  given  in 
Table  38).  At  189°C,  best  fit  is  achieved  when  EHj  is  set  at 
3*5,  giving  AAD  of  3.1%.  At  349°C,  the  optimum  EHj  is  11.0, 
with  AAD  of  2.1%.  These  experimental  data  are  correlated 
equally  well  by  the  modified  PR  equation  and  the  modified 
SRK  method,  as  illustrated  on  Figure  26.  Figure  27  compares 
the  hydrogen  K-values  calculated  by  the  modified  PR  and  the 
original  PR  equations  of  state.  At  349°C  the  experimental 
data  are  better  correlated  by  the  modified  PR  equation.  This 
is  further  evident  from  the  much  lower  AAD  values  reported 
in  Table  11  for  the  modified  PR  equation. 

In  contrast  to  the  results  of  El-Twaty  and  Prausnitz 
(14),  negative  values  of  EHj  are  found  for  some  pure 
compounds  (e.g.  quinoline)  and  most  complex  mixtures.  The 
results  for  the  K-values  of  hydrogen  in  Athabasca  Bitumen 
are  shown  on  Figure  28  for  200°  and  300°C  (the  VLE  data  are 
given  in  Table  41).  A  good  agreement  between  the 
experimental  and  the  calculated  hydrogen  K-values  is 
obtained  when  E  Hj  is  set  at  -8.6  at  200°C  and  -19.9  at 
300  °C .  The  AAD  values  are  as  lew  as  0.3%  at  200°C  and  0.4% 
at  300  0  C .  For  hydrogen-Cold  Lake  Bitumen  binary,  EHj  of  -6.6 
at  200  °C  and  -21.1  at  300°C  are  found,  giving  AAD  of  1.4% 
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FIG  28  HYDROGEN  K-VRLUE  VS  PRESSURE  FOR 


RTHfiBRSCR  BITUMEN 


HYDROGEN  K-VALUE 
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FIG  29  HYDROGEN  K-VALUE  VS  PRESSURE  FOR 


EMR  GAS  OIL 
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and  1.6%,  respectively.  The  results  for  EMR  Gas  Oil  are 
shown  on  Figure  29  (the  VLE  data  are  given  in  Table  42).  The 
optimum  EHj  of  -2.1  at  200°C  ( AAD  of  0.5%)  and  0.1  at  300°C 
( AAD  of  1.5%)  are  found  for  this  system. 

In  general ,  the  values  of  EHj  in  the  modified  PR 
equation  are  always  smaller  than  those  in  the  modified  SRK 
equation.  All  the  binary  parameters  are  shown  on  Figure  30 
where  EHj  is  plotted  against  solvent  reduced  temperature. 
The  simple  correlation  found  by  El-Twaty  and  Prausnitz  (14) 
for  seven  hydrocarbons  does  not  apply  to  most  of  the  systems 
studied  here. 

5.1.3  Grayson-Streed  Method 

The  Grayson-Streed  correlation  was  also  applied  to  the 
bicyclohexyl  (56),  tetralin  (61),  diphenylmethane  (62), 

1  -methylnaphthalene  (72),  /77-cresol  (64),  Athabasca  Bitumen 
and  Cold  Lake  Bitumen  systems.  The  results  for  the  K-values 
of  hydrogen  for  the  hydrogen-bicyclohexyl  system  are  shown 
on  Figure  31  (the  VLE  data  are  given  in  Table  37).  Despite  a 
significant  improvement  in  the  correlation  over  the 
Chao-Seader  method,  the  K-values  are  still  consistently 
underpredicted.  At  349°C,  the  AAD  is  as  high  as  24%.  For  the 
hydrogen- 1 -methylnaphthalene  system,  the  average  deviation 
is  41%  at  349 0 C . 

The  consequence  of  an  important  assumption  made  in  the 
development  of  the  Chao-Seader  and  Grayson-Streed 
correlations  explains,  in  part,  the  poor  performance  of  GS 
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procedure  for  hydrogen-heavy  hydrocarbon  systems  at  elevated 
temperatures.  Two  distinct  models  are  utilised  in  the 
development  of  the  GS  procedures;  the  'regular'  solution 
model  is  assumed  for  liquid  solution,  and  the  Redl ich-Kwong 
equation  of  state  is  used  to  predict  the  properties  of  the 
vapor  phase.  At  temperatures  approaching  the  critical  values 
the  two  phases  approach  a  dense  single  phase.  The  two  models 
will  not  necessarily  match  the  properties  for  this  phase. 

5.1.4  Correlation  of  C02  and  H2S  Systems 

The  solubility  data  for  carbon  dioxide  in  tetralin 
(57),  Athabasca  Bitumen  and  EMR  Gas  Oil  were  successfuly 
correlated  by  the  PR  equation  for  temperatures  ranging  from 
50°  to  392 0  C ,  and  pressures  up  to  6.4  MPa.  Unlike  hydrogen 
systems,  the  phase  compositions  were  strong  functions  of  the 
binary  interaction  coefficients,  6  j  j .  For  these  systems  it 
was  possible  to  find  a  value  of  the  interaction  coefficient 
that  reproduced  the  experimental  data  well  and  lay  in  the 
range  -1  < 6 ;  j  <  +1.  For  the  carbon  diox ide-tet ral in  data  of 
Sebastian  et  al.  (57)  the  interaction  coefficients  were  0.14 
at  189°  and  270°C,  0.19  at  350°C,  and  0.26  at  392°C.  The 
results  for  the  K-values  of  carbon  dioxide  are  shown  on 
Figure  32.  The  VLE  data  are  given  in  Table  43  of  Appendix 
II . 

For  the  solubility  of  carbon  dioxide  in  Athabasca 
Bitumen,  the  optimum  binary  interaction  coefficients  were 
found  to  be  0.09  at  50°,  0.07  at  100°  and  -0.04  at  300°C. 
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The  AAD  values  were  4.8%,  2.9%  and  1.1%,  respectively.  The 
results  for  the  K-values  of  carbon  dioxide  at  50°  and  300°C 
are  shown  on  Figure  33.  The  VLE  data  are  listed  in  Table  44 
of  Appendix  II.  For  the  EMR  Gas  Oil  system  the  interaction 
coefficients  were  fairly  constant  at  0.13,  0.12  and  0.13  at 
50°,  100°  and  150°C,  respectively.  The  results  for  the 
K-values  of  carbon  dioxide  are  shown  in  Figure  34.  The  data 
are  given  in  Table  45  of  Appendix  II. 

Reliable  and  accurate  phase  equilibrium  predictions  for 
hydrogen  sulphide  systems  were  also  possible  by  the  PR 
equation.  For  the  hydrogen  sulphide-Athabasca  Bitumen 
system,  optimum  interaction  coefficients  were  -0.05  at  200°C 
and  -0.055  at  300°C,  with  average  deviation  of  3.0%  and 
2.5%,  respectively.  The  dependence  of  the  K-values  of 
hydrogen  sulphide  on  pressure  are  shown  on  Figure  35.  The 
data  are  given  in  Table  46  of  Appendix  II.  For  the  EMR  Gas 
Oil  system,  the  interaction  parameters  were  0.05  at  100°C, 
and  0.18  at  200°  and  300°C.  The  results  for  the  K-values  of 
hydrogen  sulphide  are  shown  on  Figure  36.  The  data  are  given 
in  Table  47  of  Appendix  II. 

A  reasonably  good  fit  of  the  hydrogen  sulphide-n-decane 
data  of  Reamer  et  al.  (48)  was  possible  by  assigning  a 
constant  value  of  0.04  for  the  interaction  coefficient  at 
4°,  38°,  71°,  104°,  138°  and  171°C,  with  an  overall  AAD  of 
4.3%.  For  the  n~ nonane  data  of  Eakin  and  DeVaney  (12) 
overall  interaction  coefficient  of  0.04  was  also  required 
for  data  at  38°,  93°  and  204°C. 
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FIG  32  CARBON  DIOXIDE  K-VflLUE  VS  PRESSURE 


FOR  TETRALIN  (SEBASTIAN  ET  AL( 57 )  ) 


CARBON  DIOXIDE  K-VALUE 
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PRESSURE  .  MPR 


FIG  33  CBRBON  DIOXIDE  K-VBLUE  VS 


t - r 


j _ t 
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PRESSURE  FOR  RTHRBRSCR  BITUMEN 


CARBON  DIOXIDE  K-VALUE 


102 


PRESSURE  ,  MPR 

FIG  34  CRRBON  DIOXIDE  K-VRLUE  VS 


PRESSURE  FOR  EMR  GRS  OIL 


HYDROGEN  SULPHIDE  K-VflLUE 
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PRESSURE  ,  MPR 

FIG  35  HYDROGEN  SULPHIDE  K-VRLUE  VS 
PRESSURE  FOR  RTHRBRSCR  BITUMEN 


HYDROGEN  SULPHIDE  K-VHLUE 
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FIG  36  HYDROGEN  SULPHIDE  K-VflLUE  VS 


PRESSURE  FOR  EMR  GRS  OIL 
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It  is  apparent  from  the  above  results  that  the  original 
PR  equation  is  capable  of  handling  carbon  dioxide  and 
hydrogen  sulphide  in  a  wide  range  of  solvents  and  over  a 
wide  range  of  temperatures  and  pressures.  The  binary 
interaction  coefficients  are  low  and  positive,  with  the 
exception  for  carbon  dioxide  in  Athabasca  Bitumen  at  300°C 
and  hydrogen  sulphide  in  Athabasca  Bitumen  at  200°  and 
300  0 C . 

5.1.5  Fraction  Characterization 

Table  1  compares  the  pure  component  properties  required 
for  the  PR,  SRK  and  the  Grayson-St reed  methods.  It  can  be 
seen  from  this  table  that  the  PR  and  SRK  (and  their 
modifications  discussed  in  this  work)  methods  require  fewer 
pure  component  properties  than  the  Grayson-St reed 
correlation.  However,  the  PR  and  SRK  both  require 
interaction  coefficients  for  each  binary  pair  present  in  the 
system. 

For  well-defined  components  the  literature  values  of 
the  parameters  are  utilised.  In  the  cases  of  EMR  Gas  Oil  and 
the  bitumens,  the  method  of  Cavett  (3)  is  used  .  This 
requires  specification  of  normal  boiling  point  and  specific 
gravity  for  the  component. 

The  boiling  point  curve  for  the  EMR  Gas  Oil  was 
obtained  by  simulated  distillation  and  checked  against  the 
spinning-band  distillation  method.  The  molecular  weight  of 
this  solvent  was  determined  by  the  method  outlined  in  the 
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API  Data  Handbook.  The  results 

listed  below. 

for  the 

EMR  Gas 

Molecular  Weight 

250 

Specific  Gravity 

0.915 

Critical  Temperature, 

K 

793 

Critical  Pressure,  psia 

263.9 

Acentric  Factor 

0.630 

Solubility  Parameter 

( cal/cc )  0>5 

7.98 

A  boiling  point  curve  was  not  available  for  the 
bitumens  and  consequently  characterization  of  these  heavy 
hydrocarbons  was  difficult.  Molecular  weights  of  500  for 
Athabasca  Bitumen  and  475  for  Cold  Lake  Bitumen  were 
estimated  from  the  works  of  Selucky  et  al.  (58).  Specific 
gravities  were  1.002  and  1.00  for  Athabasca  and  Cold  Lake 
Bitumen,  respectively.  Cavett's  method  was  then  used  to 
derive  the  other  parameters,  assuming  a  boiling  point  of 
500°C  for  Athabasca  Bitumen  and  475°C  for  Cold  Lake  Bitumen. 
The  results  are  listed  below. 


Athabasca  Cold  Lake 


Molecular 

■  Weight 

500 

475 

Spec i f ic 

Gravity 

1  .002 

1.00 

Critical 

Temperature,  K 

947 

929 

Critical 

Pressure,  psia 

158.5 

171.7 

Acentric 

Factor 

0.963 

0.893 

' 
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Solubility  Parameter  (cal/cc)  0-5  7.365  7.3  18 


The  critical  constants  and  the  acentric  factors  of 

hydrogen ,  hydrogen  sulphide  and  carbon  dioxide  are  listed 
below . 


h2 

H2S 

co2 

Cr i t ical 

Pressure,  psia 

305 

1297 

1070 

Critical 

Temperature,  K 

41.67 

373.2 

304.2 

Acentric 

Factor 

0.0 

0.100 

0.225 

The  properties  of  hydrogen  sulphide  and  carbon  dioxide  are 
obtained  from  Reid  et  al.  (50).  The  effective  critical 
properties  of  hydrogen  are  obtained  from  Graboski  and 
Daubert  (17). 


. 


6.  Summary  and  Conclusions 

1.  A  batch  -  type  autoclave  was  installed  and  tested  for 
measurement  of  hydrogen  solubility  in  heavy  oils  and 
coal  ~  bitumen  slurries.  Solubility  data  for  hydrogen  in 
bitumens,  cracked  fractions,  heavy  gas  oil  and  coal  - 
bitumen  slurries  were  obtained  mainly'at  200°  and  300°C, 
and  pressures  up  to  24.8  MPa.  Limited  carbon  dioxide  and 
hydrogen  sulphide  solubility  measurements  were  also 
made . 

2.  The  solubility  of  hydrogen  in  pure  solvents  increased 
linearly  with  pressure,  and  increased  with  rising 
temperatures.  The  solubility  decreased  with  increase  in 
solvent  molecular  weight. 

3.  The  solubility  in  the  10  wt%  coal  -  bitumen  slurry  was 
essentially  the  same  as  in  the  bitumen.  Increase  in 
solubility  was  observed  in  the  25  wt%  coal  -  bitumen 
slurry,  and  in  the  40  wt%  coal  -  bitumen  slurry  at  200° 
and  250°C.  The  solubility  in  the  40  wt%  coal-bitumen 
slurry  at  300°C  was  lower  than  that  in  the  bitumen. 
There  is  no  explanation  available  for  this  phenomenon. 
The  type  of  coal  did  not  have  an  appreciable  effect  on 
the  hydrogen  solubility. 

4.  Carbon  dioxide  and  hydrogen  sulphide  solubility 
decreased  with  increase  in  temperature  and  with  increase 
in  solvent  molecular  weight. 

5.  The  prediction  and  correlation  of  the  solubility  data  by 
the  Peng-Robinson  equation  and  the  Gray son-St reed  method 


108 


■ 

. 

' 


109 


were  attempted.  The  use  of  the  Peng-Robinson  equation 
led  to  unrealistically  large  negative  values  of  the 
interaction  coefficient  for  some  of  the  hydrogen  - 
hydrocarbon  systems.  The  Grayson-St reed  method 
consistently  underpredicted  the  K-values  of  hydrogen. 
The  best  result  for  the  correlation  of  hydrogen 

solubility  data  was  obtained  with  the  use  of  either  a 
modified  Peng-Robinson  or  a  modified  Soave-Redl ich-Kwong 
equation  of  state,  which  introduced  an  interaction 
parameter  in  the  covolume  factor,  b. 

6.  The  Peng-Robinson  equation  was  suitable  for  the 

correlation  of  carbon  dioxide  and  hydrogen  sulphide 
solubility  data. 

7.  In  summary,  the  modified  PR  and  the  modified  SRK 

equations  of  state  work  equally  well  in  the  prediction 
of  hydrogen-hydrocarbon  VLE .  When  no  EHj  values  are 
available,  then  Figure  30  can  be  used  to  provide  an 
estimate  of  these  for  use  in  the  modified  PR  equation  of 
state.  Very  little  information  on  the  hydrogen  sulphide 
and  carbon  dioxide  systems  is  available  to  make 

estimates  of  5  j  j  for  these  gases  in  other  solvents.  The 
results  show  that  non-zero  values  have  to  be  assigned 
for  accurate  prediction  of  VLE  of  these  systems  by  the 
PR  equation  of  state. 
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TABLE  12 


1  17 


SOLUBILITY  OF  H2  IN  MESITYLENE 


T  (°  C) 

P  (kPa) 

S  (g  H2/g  Solvent) 

200 

1668 

2.74 

2082 

3.33 

3358 

5.00 

5702 

9.90 

10632 

17.73 

TABLE  13 


SOLUBILITY 


T  (°  C) 


200 


OF  H2  IN  ATHABASCA  BITUMEN  (4-AB-77) 

P  (kPa)  S  (g  H2/g  Solvent)  x104 


783 

0.574 

2024 

1  .28 

3127 

1  .93 

3955 

2.59 

5334 

3.30 

6127 

3.86 

7506 

4.78 

8884 

5.26 

9195 

4.98 

1  1607 

6.94 

13883 

9.18 

16710 

10.81 

18779 

10.90 

21158 

12.43 

24845 

15.17 

.  - 


• 

TABLE  13  continued 


1  19 


300  1187 

1.31 

2024 

2.17 

2852 

2.41 

3954 

3.80 

4402 

3.97 

5815 

4.94 

7539 

6.20 

7953 

7.29 

8469 

7.77 

1  1470 

10.08 

16710 

13.92 

1  9295 

16.93 

240  18 

20.75 

TABLE  14 


SOLUBILITY 


T  (°  C) 


200 


300 


OF  H2  IN  COLD  LAKE  BITUMEN  (1-CL-77) 

P  (kPa)  S  (g  H2/g  Solvent)  x104 


1920 

1  .25 

2037 

0.89 

2988 

1  .92 

5506 

3.36 

8402 

4.88 

1  3261 

8.13 

1402  1 

8.30 

1  7089 

10.29 

20916 

12.89 

23983 

14.77 

2300 

2.06 

4367 

4.28 

8505 

8.04 

1  2849 

12.76 

1  7468 

16.73 

22776 

18.80 

I 

■ 
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SOLUBILITY 


T  (°  C) 


200 


300 


TABLE  15 


OF  H2  IN  LLOYDMINSTER  RESIDUE  (3-LL-77) 

P  (kPa)  S  (g  H2/g  Solvent)  x104 


2058 

0.94 

3815 

1  .79 

7125 

3.67 

1  0298 

5.29 

14020 

7.64 

1  7227 

9.28 

21226 

11.88 

24603 

13.75 

1440 

0.92 

3196 

2.89 

3417 

1  .77 

5354 

3.60 

6919 

4.90 

9746 

7.67 

1  2086 

8.99 

14778 

10.81 

19295 

15.60 

23052 

18.88 

24259 

17.64 

• 

TABLE  1 6 


SOLUBILITY  OF  H2  IN  650°  F+  TOPPED 
ATHABASCA  BITUMEN  (2-AB-77) 

T  (°  C)  P  (kpa)  S  (g  H2/g  Solvent)  x104 


200 

3334 

1  .77 

8905 

4.79 

15179 

7.64 

19316 

10.30 

300 

3610 

3.40 

5808 

4.74 

8878 

6.62 

13980 

11.48 

19757 

15.08 

370 

4293 

7.05 

638  1 

12.63 

7706 

15.18 

1  1470 

23.79 

13165 

29.02 

123 


TABLE  17 

SOLUBILITY  OF  H2  IN  800°  F+  TOPPED 
ATHABASCA  BITUMEN  (36-1) 

T  (°  C)  P  (kPa)  S  (g  H2/g  Solvent)  x104 


200 

3044 

1.58 

1  0534 

5.12 

1  6294 

8.06 

22534 

11.39 

300 

5594 

3.12 

13074 

7.59 

22034 

12.41 

370 

4204 

2.79 

1  1463 

8.52 

18084 

12.70 

2251  3 

16.00 

■ 

TABLE  18 


SOLUBILITY  OF  H2  IN  650°  F+  TOPPED 
COLD  LAKE  BITUMEN  (37-1) 


(°  c) 

P  ( kPa ) 

S  (g  H2/g  Solvent) 

200 

5693 

2.94 

1  0834 

5.91 

18294 

9.55 

22894 

12.47 

300 

3504 

2.53 

10493 

8.06 

17154 

13.11 

23192 

18.37 

370 

3543 

6.90 

5483 

14.32 

8593 

25.01 

1  1564 

36.69 

*  . 


TABLE  19 


SOLUBILITY  OF  H2  IN  800°  F+  TOPPED 
COLD  LAKE  BITUMEN  (6-CL-77) 

T  (°  C)  P  (kPa)  S  (g  H2/g  Solvent)  x104 


5084 

3.26 

1  1663 

6.40 

1  7444 

9.61 

24744 

13.81 

5064 

3.88 

1  1043 

7.63 

15494 

11.05 

23403 


17.07 


TABLE  20 


SOLUBILITY  OF  H2  IN  650°  F+ 
LLOYDMI NSTER  RESIDUE  (3-LL-78) 


(°  c) 

P  ( kPa ) 

S  (g  H2/g  Solvent) 

200 

3596 

2.39 

1  1697 

7.50 

20191 

12.63 

300 

1611 

1  .46 

4437 

3.72 

9058 

7.02 

1  5263 

13.07 

2002  1 

17.12 

127 


TABLE  21 

SOLUBILITY  OF  H2  IN  760°  F+ 
LLOYDMINSTER  RESIDUE  (1-LL-78) 

T  (°  C)  P  (kPa)  S  (g  H2/g  Solvent)  x104 


200 

991 

0.83 

2575 

1  .60 

2886 

1  .83 

5099 

2.71 

5713 


3.61 


TABLE  22 


SOLUBILITY  OF  H2  IN  AMOCO  FEED  (4-AM-78) 


(°  c) 

P  (kPa) 

S  (g  H2/g  Solvent) 

200 

5614 

4.76 

1  1893 

9.87 

17753 

14.67 

23594 

19.04 

300 

5744 

6.61 

11194 

12.79 

17394 

20.00 

20594 


22.31 


' 


• 

■ 

TABLE  23 


SOLUBILITY  OF  H2  IN  EMR' GAS  OIL 


(°  c) 

P  (kPa) 

S  (g  H2/g  Solvent) 

200 

3333 

3.16 

7403 

7.03 

9452 

9.08 

12935 

12.31 

18843 

17.71 

21643 

20.15 

300 

3694 

4.60 

9844 

12.60 

14293 

17.56 

18814 

24.52 

TABLE  24 


SOLUBILITY  OF  H2  IN  HEAVY  ENDS  (91-1-1) 


T  (°  C) 

P  ( kPa ) 

S  (g  H2/g  Solv 

100 

2733 

1  .54 

7033 

3.74 

1  3743 

7.27 

21453 

10.77 

200 

3294 

2.53 

10114 

7.58 

17044 

12.59 

23544 

18.33 

23954 

18.77 

300 

3974 

3.87 

9744 

9.92 

16864 

16.68 

23264 

22.92 

TABLE  25 


SOLUBILITY  OF  H2  IN  HEAVY  ENDS  (78-T-22) 

T  (°  C)  P  (kPa)  S  (g  H2/g  Solvent)  xIO4 


3941 

5.48 

8659 

14.18 

11981 

18.82 

1  2879 

20.02 

18116 

27.35 

TABLE  26 


SOLUBILITY  OF  H2  IN  HEAVY  ENDS  (78-T-23) 

T  (°  C)  P  (kPa)  S  (g  H2/g  Solvent)  xIO4 


200 

6114 

4.35 

1  1694 

8.12 

16914 

11.33 

20594 

13.90 

300 

5402 

5.12 

1  1933 

10.78 

18193 

16.80 

23993 

21.11 

• 

TABLE  27 


SOLUBILITY  OF  H2  IN  HEAVY  ENDS  (78-T-24) 

T  (°  C)  P  (kPa)  S  (g  H2/g  Solvent)  xIO4 


200 

4993 

3.56 

12103 

8.66 

18604 

12.69 

24444 

16.29 

300 

4484 

4.60 

10134 

9.61 

16134 

15.43 

22534 


21.14 


• 
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TABLE  28 

SOLUBILITY  OF  H2  IN  COAL-BITUMEN  SLURRY 


10 

WT  %  SUB-BITUMINOUS  COAL 

90  WT  5 

l  COLD  LAKE  BITUMEN  (1-CL-77) 

T  (°  C) 

P  (kPa)  S  (g  H2/g  Solvent)  x104 

200 

4243  2.32 

8599  5.11 

12944  7.88 

18894  11.49 

21844  14.40 

24344  15.88 

300 

4844  3.86 

8043  7.09 

8683  7.56 

14833  12.36 

20303  17.79 

23182  20.86 

135 


TABLE  29 

SOLUBILITY  OF  H2  IN  COAL-BITUMEN  SLURRY 
40  WT  %  SUB- BITUMINOUS  COAL 
60  WT  %  COLD  LAKE  BITUMEN  (1-CL-77) 


T  (°  C)  P  (kPa) 


200  4284 

6154 
8888 
10044 
13143 
1  7773 
20343 
22694 

250  5233 

7445 
8972 
1  1263 
1  1553 
1  6594 
20524 


S  (g  H2/g  Solvent)  x104 


2.82 

4.22 

6.37 

7.68 

9.73 

13.32 

14.78 

16.57 

1  .80 
3.45 
4.86 
7.26 
8.20 
12.12 
16.03 


TABLE  29  continued 


300  4494  0.603 

6844  0.995 

6695  1.14 

8393  2.41 

8543  2.35 

10134  3.31 

11404  4.75 

12292  5.88 

14242  8.71 

14452  8.79 

15843  11.12 

16144  12.28 

17623  13.67 

19044  16.79 

20513  17.63 

21993  20.01 


TABLE  30 


SOLUBILITY  OF  H2  IN  COAL-BITUMEN  SLURRY 
40  WT  %  LIGNITE  COAL 
60  WT  %  COLD  LAKE  BITUMEN  (1-CL-77) 

T  (°  C)  P  ( kPa )  S  (g  H2/g  Solvent)  x104 


200 

6342 

2.74 

7993 

6.12 

12473 

9.04 

16982 

12.40 

250 

4793 

2.09 

5584 

2.79 

9154 

5.75 

13424 

10.03 

15585 

11.36 

20394 

15.75 

300 

6473 

0.739 

7544 

1.11 

9553 

3.80 

11615 

6.11 

14015 

8.44 

17593 

13.66 

17.27 


21043 


. 


TABLE  31 


SOLUBILITY  OF  H2  IN  COAL-BITUMEN  SLURRY 
25  WT  %  SUB: BITUMINOUS  COAL 
75  WT  %  COLD  LAKE  BITUMEN  (1-CL-77) 

T  (°  C)  P  (kPa)  S  (g  H2/g  Solvent)  xIO4 


200 

5543 

3.17 

8844 

5.66 

12043 

7.88 

16283 

10.92 

20492 

13.70 

300 

4144 

4.03 

7225 

7.25 

10195 

10.24 

11013 

9.25 

16653 


15.52 


. 
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TABLE  32 


SOLUBILITY  OF  H2S  IN  ATHABASCA  BITUMEN  (4-AB-77) 


T  (°  C)  P  ( kPa )  S  (g  H2S/g  Solvent)  xIO3 

200  390.4  2.65 

773.0  5.99 

1087.0  8.85 

1390.0  11.80 

1624.0  13.30 

2127.0  19.57 

300  389.7  0.99 

404.7  1.83 

624.1  1.91 

659.0  2.52 

865.7  3.87 

862.4  4.38 

1010.0  4.08 

1238.0  6.35 

1300.0  6.93 

1452.0  7.88 

1652.0  7.48 

1745.0  8.52 

1824.0  10.55 


• 

TABLE  32  continued 


2061.0 

11.10 

2068.0 

13.48 

2368.0 

11.26 

2458.0 

15.14 

3357.0 

21.70 

TABLE  33 


SOLUBILITY  OF  H2S  in  HEAVY  ENDS  (78-T-24) 

T  (°  C)  P  (kPa)  S  (g  H2S/g  Solvent)  xIO3 


200  273.5  1.20 

482.1  2.93 

831.1  4.77 

1224.3  10.52 

1741.9  14.48 

2534.9  24.78 

3468.6  34.63 

300  384.3  1.46 

551.7  1.46 

797.2  2.02 

1180.1  3.79 

1480.1  5.11 

1952.7  8.60 

2582.9  14.64 


* 
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TABLE  34 


SOLUBILITY  OF  H2S  IN  EMR  GAS  OIL  ' 

T  (°  C)  P  (kPa)  S  (g  H2S/g  Solvent)  xIO2 


100 

393.4 

0.823 

750.4 

1  .69 

768.4 

1  .67 

1273 

2.74 

1613 

4.33 

200 

425.3 

0.241 

780.4 

0.502 

1094 

0.811 

1469 

1  .30 

1982 

2.04 

300 

735.5 

0.274 

946.9 

0.363 

1400 

0.755 

1718 

1.10 

2136 

1  .54 

143 


TABLE  35 

SOLUBILITY  OF  C02  IN  ATHABASCA  BITUMEN  (4-AB-77) 


T  (°  C)  P  (kPa)  S  (g  C02/g  Solvent)  xIO2 

50  1079  1.17 

1479  1.92 

2279  2.82 

2962  3.58 

3788  4.42 

4461  5.66 

5453  6.00 

100  1535  1.16 

3037  2.28 

3810  2.88 

4568  3.58 

5471  4.90 

300  1858  0.750 

2859  1.17 

4279  1.83 

5689  2.64 

6444  2.92 


TABLE  36 


SOLUBILITY  OF  C02  IN  EMR  GAS  OIL 


T  (  0  C) 

P  (kPa) 

S  (g  C02/g  Solvent) 

50 

635.0 

0.853 

707.2 

1.14 

1589 

2.70 

301  1 

5.30 

4362 

8.42 

5637 

12.46 

100 

845.6 

0.843 

2221 

2.47 

3693 

4.42 

483  1 

5.76 

5596 

6.97 

150 

755.9 

0.512 

1604 

1  .22 

2399 

1  .98 

3540 

2.89 

5084 

5.22 

6188 

6.90 

. 

• 
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TABLE  37 


VAPOR  -  LIQUID  EQUILIBRIUM  DATA 
FOR  H2  -  BI CYCLOHEXYL  (56) 

Liq.  Comp,  of  H2  Equil.  Const,  of  H2 

P(kPa)  Exp  Mod  SRK  GS  Exp  Mod  SRK  GS 


189° 

C ( E  H j  =7 . 

4 (Mod  : 

SRK)  ) 

2027 

0.0181 

0.0174 

0.0201 

54.4 

56.6 

48.9 

304  1 

0.0262 

0.0240 

0.0257 

37.7 

38.  1 

32.9 

5068 

0.0443 

0.0425 

0.0496 

22.4 

20.2 

23.4 

10135 

0.1158 

0.1157 

0.1389 

8.61 

8.62 

7.18 

20270 

0.1418 

0.1473 

0.1786 

7.03 

6.77 

5.59 

25338 

0.1718 

0.1761 

0.2153 

5.81 

5.67 

4.64 

348° 

C ( E  H  j  =2 1 

. 4 (Mod 

SRK)  ) 

2027 

0.0211 

0.0227 

0.0230 

29.0 

26.0 

25.6 

304  1 

0.0382 

0.0393 

0.0404 

19.16 

18.17 

17.61 

5068 

0.0686 

0.0644 

0.0691 

12.09 

11.60 

10.95 

10135 

0.1398 

0.1399 

0.1559 

6.48 

6.45 

5.77 

1  5203 

0.1975 

0.1986 

0.2327 

4.71 

4.69 

4.00 

20270 

0.2578 

0.2490 

0.3056 

3.66 

3.79 

3.11 

25338 

0.3112 

0.2929 

0.3752 

3.06 

3.26 

2.57 

. 


■ 

. 

TABLE  38 


VAPOR  -  LIQUID  EQUILIBRIUM  DATA 
FOR  H2  -  B I CYCLOHEXYL  (56) 


Lig. 

Comp,  of 

h2 

Equi 1 . 

Const 

.  of  H 

P(kPa) 

Exp 

PR 

Mod  PR 

Exp 

PR 

Mod  PR 

189°  C ( 6 i 

j=-1 . 4 ( PR ) f  E  H j  =3 

. 5 (Mod  PR) ) 

2027 

0.0181 

0.0172 

0.0172 

54.4 

57.2 

57.3 

304  1 

0.0262 

0.0257 

0.0257 

37.7 

38.5 

38.5 

5068 

0.0442 

0.0420 

0.0421 

22.4 

23.6 

23.6 

10135 

0.0810 

0.0804 

0.0808 

12.29 

12.38 

12.32 

15203 

0.1158 

0.1153 

0.1162 

8.61 

8.65 

8.58 

20270 

0.1418 

0.1470 

0. 1488 

7.03 

6.78 

6.70 

25338 

0.1718 

0.1762 

0.1788 

5.81 

5.66 

5.58 

348°  C ( 6  j 

j=-1 . 4 ( PR ) , E  H j  =  1 1 

. 0 (Mod 

PR)  ) 

2027 

0.0211 

0.0207 

0.0221 

29.0 

28.9 

26.7 

304  1 

0.0382 

0.0356 

0.0383 

19.16 

20.3 

18.67 

5068 

0.0688 

0.0635 

0.0691 

12.09 

13.07 

11.86 

10135 

0.1398 

0.1240 

0.1382 

6.48 

7.38 

6.52 

15203 

0.1975 

0.1740 

0.1978 

4.71 

5.42 

4.70 

20270 

0.2578 

0.2159 

0.2499 

3.66 

4.44 

3.77 

25338 

0.3112 

0.2517 

0.2959 

3.06 

3.84 

3.21 

, 


• 

TABLE  39 


VAPOR  -  LIQUID  EQUILIBRIUM  DATA 
FOR  H2  -  TETRALIN  (61 ) 


Liq. 

Comp.  ' 

of  H2 

Equi 1 . 

Const 

.  of  H 

P(kPa)  Exp 

PR 

Mod  PR 

Exp 

PR  ] 

Mod  PR 

169.6°  C 

(5, ,=-0 

. 6 ( PR) , E  H j 

=4 . 6 (Mod 

PR)  ) 

2027 

0.0118 

0.0119 

0.0117 

81.8 

81.0 

82.4 

304  1 

0.0176 

0.0179 

0.0176 

55.4 

54.5 

55.3 

5068 

0.0297 

0.0296 

0.0292 

33.  1 

33.3 

33.7 

10135 

0.0571 

0.0570 

0.0567 

17.33 

17.40 

17.48 

15203 

0.0823 

0.0822 

0.0823 

12.06 

12.09 

12.07 

20270 

0.1051 

0.1055 

0.1062 

9.46 

9.44 

9.37 

25338 

0.1289 

0.1270 

0.1285 

7.72 

7.85 

7.74 

348.6°  C 

(6 ,  j  =-0 

. 6 ( PR ) , E  H j 

=13. 8 (Mod 

PR)  ) 

5068 

0.0452 

0.0487 

0.0452 

15.63 

14.66 

15.40 

10135 

0.0925 

0.1019 

0.0975 

9.02 

8.36 

8.53 

1  5203 

0.1390 

0.1467 

0.1436 

6.35 

6.14 

6.12 

20270 

0.1884 

0.  1849 

0.1847 

4.80 

5.00 

4.89 

25338 

0.2314 

0.2179 

0.2216 

3.96 

4.31 

4.  14 

• 

. 

K*  <:.r  . 

TABLE  40 


VAPOR  -  LIQUID  EQUILIBRIUM  DATA 
FOR  H2  -  TETRALIN  (61 ) 


Liq. 

Comp,  of 

h2 

Equil 

.  Const. 

of  H 

P(kPa) 

Exp 

Mod  SRK 

GS 

Exp 

Mod  SRK 

GS 

169. 

6°  C ( E  H j  = 

7. 98 (Mod 

SRK)  ) 

2027 

0.0118 

0.0118 

0.0145 

81.8 

81.9 

66.7 

304  1 

0.0176 

0.0177 

0.0218 

55.4 

55.  1 

44.8 

5068 

0.0297 

0.0292 

0.0362 

33.  1 

33.7 

27.2 

10135 

0.0571 

0.0564 

0.0705 

17.33 

17.59 

14.05 

1  5203 

0.0823 

0.0813 

0.1027 

12.06 

12.23 

9.68 

20270 

0.1051 

0.1042 

0.  1326 

9.46 

9.55 

7.50 

25338 

0.  1289 

0.1254 

0.1605 

7.72 

7.94 

6.21 

348. 

6°  C  ( E  Hj 

=22.4 (Mod 

SRK)  ) 

5068 

0.0452 

0.0465 

0.0524 

15.63 

14.97 

13.21 

10135 

0.0925 

0.0991 

0.1176 

9.02 

8.41 

6.83 

15203 

0.1390 

0.  1446 

0.1800 

6.35 

6.10 

4.86 

20270 

0.1884 

0.1841 

0.2403 

4.80 

4.92 

3.76 

25338 

0.2314 

0.2188 

0.2988 

3.96 

4.21 

3.09 

• 

TABLE  41 


VAPOR  -  LIQUID  EQUILIBRIUM  DATA  FOR 
H2  -  ATHABASCA  BITUMEN  (4-AB-77) 

Liq.  Comp,  of  H2  Equil.  Const  of  H 


P(kPa) 

Exp 

Mod  PR 

Exp 

Mod  PR 

200°  C 

(EHj=-8.6) 

• 

2024 

0.0307 

0.0302 

32.6 

33.  1 

5334 

0.0757 

0.0755 

13.21 

13.24 

1  1607 

0.1468 

0.1497 

6.8  1 

6.68 

16710 

0.2111 

0.2010 

4.74 

4.98 

21158 

0.2356 

0.2405 

4.24 

4.16 

24845 

0.2734 

0.2701 

3.66 

3.70 

300°  C 

(E  Hj  =-  1  9 . 9) 

2852 

0.0564 

0.0597 

17.73 

16.73 

4402 

0.0896 

0.0889 

11.16 

11.25 

7953 

0.1531 

0. 1484 

6.53 

6.74 

1  1470 

0.2000 

0.1993 

5.00 

5.02 

16710 

0.2566 

0.2638 

3.90 

3.79 

240  18 

0.3398 

0.3369 

2.94 

2.97 

. 


' 

■ 

TABLE  42 


VAPOR 

-  LIQUID 

'  EQUILIBRIUM 

DATA 

FOR  H2  - 

EMR  GAS  OIL 

Liq.  Comp,  of  H2 

Equil.  < 

Const  of  H 

P(kPa) 

Exp. 

Mod  PR 

Exp. 

Mod  PR 

200° 

C  ( E  H j  =~3 . 2 ) 

3333 

0.0377 

0.0378 

26.5 

26.4 

7403 

0.0802 

0.0801 

12.47 

12.47 

9452 

0.1012 

0.0999 

9.88 

10.00 

12935 

0.1325 

0.1316 

7.55 

7.59 

18843 

0.1801 

0.1803 

5.55 

5.54 

21643 

0.1999 

0.2014 

5.00 

4.96 

300° 

C  ( E  H  j  =  0  .  1  ) 

3694 

0.0540 

0.0543 

18.52 

18.04 

9844 

0.1351 

0.1340 

7.40 

7.40 

14293 

0.1788 

0.1837 

5.59 

5.41 

18814 

0.2332 

0.2286 

4.29 

4.35 

. 

TABLE  43 


VAPOR  -  LIQUID  EQUILIBRIUM  DATA 
C02  -  TETRALIN  (57) 


P(kPa) 

Liq.  Comp. 

Exp 

of  C02 

PR 

Equil.  Const 

Exp 

.  of  CO 

PR 

189° 

C 

(6 , i=0. 14) 

2047 

0.0555 

0.0556 

17.30 

17.21 

3085 

0.0823 

0.0836 

11.79 

11.57 

4095 

0.1095 

0.1100 

8.90 

8.84 

5128 

0.1370 

0.1364 

7.13 

7.15 

270° 

C 

(6 , j=0. 14) 

2017 

0.0429 

0.0421 

18.47 

18.45 

3055 

0.0673 

0.0679 

12.64 

12.36 

4095 

0.0931 

0.0933 

9.44 

9.32 

5098 

0.1149 

0.1172 

7.77 

7.56 

350° 

C 

(6 , j=0. 19) 

1999 

0.0229 

0.0212 

15.22 

14.97 

3069 

0.0492 

0.0494 

10.61 

10.21 

4105 

0.0764 

0.0762 

8.00 

7.82 

5057 

0.1018 

0.1009 

6.49 

6.43 

. 
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TABLE  43  Continued... 


392°  C  ( 6  j  j  =0 . 26 ) 


3118 

0.0340 

0.0338 

7.85 

7.30 

4145 

0.1174 

0.1168 

5.93 

5.75 

5108 

0.0902 

0.0932 

4.99 

4.79 

1079 

2279 

2962 

3788 

4461 

5453 

1535 

3037 

3810 

4568 

547  1 
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TABLE  44 


VAPOR 

-  LIQUID  EQUILIBRIUM  DATA 

FOR 

CO 2  -  ATHABASCA 

BITUMEN 

Liq.  Comp. 

of  CO 2 

Equil.  Const 

.  of 

Exp 

PR 

Exp 

PR 

50°  C  (6  |  j 

=0.09) 

0.1174 

0.1138 

8.52 

8.79 

0.2427 

0.2256 

4.12 

4.43 

0.2892 

0.2828 

3.46 

3.54 

0.3343 

0.3461 

2.99 

2.89 

0.3914 

0.3932 

2.48 

2.54 

0.4054 

0.4556 

2.47 

2.20 

100°  C  (5  j 

,=0.07) 

0.1165 

0.1093 

8.58 

9.15 

0.2058 

0.2039 

4.86 

4.91 

0.2466 

0.2481 

4.06 

4.03 

0.2892 

0.2887 

3.46 

3.46 

0.3577 

0.3339 

2.80 

3.00 

i.r  ‘ 

■ 

TABLE  44  Continued 


300°  C  (6  j  j  =  —  0 . 0 4 ) 


1858 

0.0785 

0.0787 

12.73 

12.70 

2859 

0.1174 

0.1185 

8.52 

8.43 

4279 

0. 1722 

0.1721 

5.81 

5.81 

5689 

0.2308 

0.2222 

4.61 

4.50 

6444 

0.2491 

0.2479 

4.14 

4.03 

TABLE  45 


VAPOR 

-  LIQUID 

EQUILIBRIUM  FOR 

CO 2  “  EMR  GAS  OIL 

Liq.  Comp 

'.  of  CO  2 

Equil.  Const 

.  of 

P(kPa) 

Exp 

PR 

Exp 

PR 

50°  C 

<  6  ,  j  =0 .  1  3 ) 

635.0 

0.0497 

0.0540 

20.  1 

18.53 

707.2 

0.0608 

0.0599 

16.44 

16.68 

1589 

0.1329 

0.1297 

7.52 

7.71 

301  1 

0.2315 

0.2313 

4.32 

4.32 

4362 

0.3237 

0.3163 

3.09 

3.16 

5637 

0.4146 

0.3865 

2.41 

2.59 

100°  C 

( 6 , j  =0 .  12) 

845.6 

0.0456 

0.0488 

21.7 

20.5 

222  1 

0.1226 

0.  1225 

8.07 

8.16 

3693 

0.1986 

0.  1942 

4.97 

5.10 

483  1 

0.2445 

0.2448 

4.07 

4.08 

5596 

0.2808 

0.2767 

3.54 

3.61 

. 


. 

• 

. 
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TABLE  45  Continued... 


150°  C  (6 ,  j=0.  13) 


755.9 

0.0282 

0.0326 

35.2 

30.7 

1604 

0.0649 

0.0676 

15.35 

14.78 

2399 

0.1005 

0.0991 

9.82 

10.08 

3540 

0.  1400 

0.1422 

7.10 

7.03 

5084 

0.2271 

0.1966 

4.36 

5.09 

6188 

0.2790 

0.2329 

3.47 

4.29 

• 

TABLE  46 


VAPOR  -  LIQUID  EQUILIBRIUM  DATA  FOR 
H2S  -  ATHABASCA  BITUMEN  (4-AB-77) 


Liq.  Comp. 

of  H2S 

Equi 1 . 

Const,  of 

P(kPa) 

Exp 

PR 

Exp 

PR 

200° 

C 

(6  j  j  =-0 . 05 ) 

585.7 

0.0556 

0.0585 

17.99 

17.10 

822.9 

0.0811 

0.0815 

12.33 

12.27 

1049 

0.1053 

0.1031 

9.50 

9.70 

1263 

0.1282 

0.1230 

7.80 

8.13 

1746 

0.  1807 

0.1667 

5.53 

6.00 

2155 

0.2273 

0.2029 

4.40 

4.93 

300° 

c 

(  6  j  j=  —  0.055) 

429.9 

0.0286 

0.0302 

35.0 

33.  1 

828.8 

0.0556 

0.0575 

17.99 

17.38 

1202 

0.0811 

0.0825 

12.33 

12.11 

1551 

0.1053 

0.1054 

9.50 

9.48 

1874 

0.1282 

0.1262 

7.80 

7.92 

257  1 

0.1807 

0.1696 

5.53 

5.89 

3109 

0.2273 

0.2020 

4.40 

4.95 

. 


■ 

TABLE  47 


VAPOR  -  LIQUID  EQUILIBRIUM  DATA 
FOR  H2S  -  EMR  GAS  OIL 


Liq.  Comp. 

of  H2S 

Equi 1 . 

Const,  of  H 

P(kPa) 

Exp 

PR 

Exp 

PR 

100°  C 

(6,  j  =0 . 05 ) 

393.4 

0.0572 

0.0579 

17.27 

17.48 

750.4 

0.1101 

0.1086 

9.21 

9.08 

768.4 

0.1090 

0.1111 

9.00 

9.17 

1273 

0.1678 

0.  1800 

5.55 

5.96 

1613 

0.2414 

0.2247 

4.45 

4.14 

200°  C 

(6,  j  =0 . 18  > 

425.2 

0.0174 

0.0233 

44.5 

57.5 

780.4 

0.0356 

0.0406 

24.5 

28.  1 

1094 

0.0563 

0.0566 

17.60 

17.76 

1469 

0.0869 

0.0753 

11.51 

13.23 

1982 

0.1304 

0.1004 

7.67 

9.93 

300°  C 

(8  i  j-0. 18) 

946.9 

0.0260 

0.0399 

38.5 

23.  1 

1397 

0.0526 

0.0597 

19.01 

15.79 

1718 

0.0748 

0.0736 

13.37 

12.92 

2136 

0.1016 

0.0915 

9.84 

10.48 

* 

• 

10.  Appendix  III 
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TABLE  48 

Hempel  Distillation  of  Athabasca  Topped  Bitumen  (2-AB-77) 


Fraction 

Cut  0 

C  Vol  % 

1 

50 

2 

75 

3 

100 

4 

125 

5 

150 

6 

175 

0.3 

7 

200 

0.4 

8 

225 

0.6 

9 

250 

1  .2 

10 

275 

6.2 

Distillation  at 

40  mm  Hg 

1  1 

200 

1  .7 

12 

225 

1  .  1 

13 

250 

3.9 

14 

275 

7.5 

15 

300 

12.2 

Residuum 

64 . 6 

Loss 

0.3 

Initial 

Boiling 

Point  152 

Sum  %  Sp.  Gravity  15/15°  C 


0.3 

0.824 

0.7 

0.824 

1  .3 

0.824 

2.5 

0.838 

8.7 

0.868 

10.4 

0.881 

11.5 

0.902 

15.4 

0.927 

22.9 

0.940 

35.  1 

1  .058 

99.7 

1  .058 

C 


: ;  . 

• 
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TABLE  49 

Hempel  Distillation  of  Lloydminster  Residue  (3-LL-78) 


Fraction 

Cut  0  C 

Vol  % 

1 

50 

2 

75 

3 

100 

4 

125 

5 

150 

6 

175 

7 

200 

8 

225 

0.3 

9 

250 

1  .2 

275 

3.0 

Distillation  at  40 

mm  Hg 

1  1 

200 

1  .2 

12 

225 

3.9 

13 

250 

6.3 

14 

275 

5.7 

15 

300 

9.2 

Residuum 

% 

69.3 

Gain 

0.  1 

Initiial 

Boiling  Point  218 

Sum  %  Sp.  Gravity  15/15°  C 


0.3 

1.5 

4.5 

5.7 

9.6 
15.9 
21.6 
30.8 

100.1 


0  C 


- 
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TABLE  50 

Hempel  Distillation  of  Amoco  Feed 


Fraction  Cut  0  C  Vol  %  Sum  %  Sp.  Gravity  15/15°  C 


1 

50 

2 

75 

3 

100 

4 

125 

5 

150 

6 

175 

7 

200 

3.3 

3.3 

0.800 

8 

225 

8.0 

11.3 

0.820 

9 

250 

9.9 

21.2 

0.843 

10 

275 

10.4 

31.6 

0.859 

Distillation  at 

40  mm  Hg 

1  1 

200 

7.7 

39.3 

0.874 

12 

225 

7.2 

46.5 

0.889 

13 

250 

3.9 

50.4 

0.912 

14 

275 

3.7 

54.  1 

0.935 

15 

300 

6.5 

60.6 

0.964 

Residuum 

% 

38.  1 

98.7 

Loss 

1  .3 

Initial 

Boiling 

Point  176° 

C 

■ 


• 

* 
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